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REMARKS 



Reconsideration of this patent application is respectfully 
requested in view of the foregoing amendments, and the following 
remarks . 

The amendments to this patent application are as follows. 
Amendments have been made to the drawings, the Specification, the 
Abstract, and the claims so as to overcome all of the various 
formal rejections by the Patent Examiner. 

The amendments to the drawings are to cancel the drawings 
now on file and to replace them with a revised set of drawings 
for FIG. 1 to 11 . These revised drawings have deleted any 
reference to a total of only 10 drawing sheets and have been 
amended to refer to a total of 11 drawing sheets. Also these 
revised drawings have used "Prior Art" to identify each one of 
FIGS. 3a, 3b, 10a and 10b. 

The Substitute Specification has been amended as follows. 
On Page 1, the Substitute Specification has been amended to 
delete any reference to the PCT patent application, as required 
by the Patent Examiner. The Substitute Specification on Page 3 0 
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has been amended in order to correct the spelling of the word 
"photometric". The Substitute Specification on Page 32 has been 
amended to refer to "figure 1" in line 6. The Substitute 
Specification on Page 39 has been amended to revise the equations 
to show larger print for the subscripts. 

The Abstract of the Disclosure has been amended to delete 
"legal phraseology" and "purported merits", objected to by the 
Patent Examiner. Thus the word "said" has been changed to "the". 
Also the language "which are useful" and "that can be used", have 
been cancelled. 

The amendments to the claims are as follows. 

The Patent Examiner had objected to claim 92 in line 3, 
because the terminology "the infrared active components" lacks 
clear antecedent basis. Thus claim 92 has now been amended to 
provide proper antecedent basis. 

Claim 96 has been amended to correct a minor typographical 
error concerning the density. 
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The Patent Examiner had objected to claim 98, in lines 17- 
18, because the terminology ''the radiation receiver'' has a 
different antecedent basis which is ''at least one radiation 
receiver'' . Thus claim 98 has now been amended to provide the 
proper antecedent basis. 

The Patent Examiner objected to claim 105, in lines 5 and 1 , 
because reference numerals were not deleted. Thus claim 105 has 
been amended in order to delete these reference numerals. 

The Patent Examiner objected to claim 108 , in lines 2-4 of 
this claim which now recites "..wherein the probe cell is 
sweepable with an infrared inactive gas, preferably for carrying 
out a null measurement..." which is still allegedly indefinite 
because it is not clear as to the range of the probe cell being 
sweeped with an infrared inactive gas. Thus claim 108 has been 
amended to delete the "preferably" terminology. 

The Patent Examiner has repeated his objection with respect 
to the terminology of industry standards being employed in claims 
87, 88, 89 and 94. 
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with regard to the formal rejection by the Patent Examiner 
regarding these claims, relating to the use of standard 
calculation methods, the following should be noted: 

The standard calculation methods used in the present patent 
application are both shown in ISO 12213, specifically AGA8-92DC 
in ISO 12213-2 and GERG in ISO 12213-3. No English version of 
the DVGW working paper is known, but it is derived from 12213. 
This means that the English language document ISO 12213 is the 
relevant document. Enclosed is a copy of ISO 12213 in English. 

At present, only the two methods mentioned are of practical 
significance, since they produce the best results. An older 
method is the method according to AGA NX19, which is allowed to 
be used only at low pressures, for example. However, a new GERG2 
method is currently being developed, which is superior to the 
existing equations and is also supposed to be elevated to a 
standard. Moreover, the time when this will happen is not known 
as yet. However, this method will differ from the previous 
method also in terms of the designation. Hence the uncertainty 
that exists from the point of view of the Patent Examiner, with 
regard to the scope of the standard referred to, in each 
instance, should actually no longer exist. It is pointed out to 
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the Patent Examiner, with regard to the enclosure, that it is 
true that new standards can be developed. These new standards 
are also being used in place of the old standards, so that even 
in the relevant art, no lack of clarity can occur with regard to 
the standard that is actually meant. For these reasons, the 
formal ground rejection with regard to the precise content of the 
claims is respectfully submitted to have been overcome. 

A restriction of the claims only to the methods mentioned is 
not necessary, since there is no prior art rejection. Also one 
skilled in the art would clearly understand the meaning of all 
the art recognized terminology. 

However, the standard calculation method, as an iteration 
method according to the present invention, is also completely 
independent of the real gas calculation method actually being 
used. Instead, any algorithm for describing the real gas 
behavior of natural gas that yields one of the following variable 
is suitable; 

real gas factor Z = compression factor Z 

compressibility K - compressibility ratio K 
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Based upon all of the above-discussed amendments to the 
Abstract, the Specification, the drawings, and all the claims , 
the present patent application is now in complete compliance with 
all the requirements of 35 U.S.C. 112. Withdrawal of this ground 
of rejection is respectfully requested. 



A prompt notification of allowability is respectfully 



requested. 



Respectfully submitted. 



COLLARD Sc ROE, P.C. 
1077 Northern Boulevard 
Roslyn, New York 11576 
(516) 365-9802 



DIETER STIRNBERG ET AL-1 PCT 



By: 



Allison C.Collard,Reg^^No.22, 532 
Edward R. Freedman, Reg .No . 26 , 048 
Attorneys for Applicant 
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Enclosures : 



1) Drawing Replacement Sheets for FIG. 1 to 1 1 . 

2) Copy of ISO 1 22 1 3 in English 



I hereby certify that this correspondence is being deposited with the U.S. Postal Service as first class mail in an 
envelope addressed to: Commissioner of Patents, P.O. Box 1450>.lexaTtfiria, VA 22313-1450, on June 23, 2004 

Maria Guastella 
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Foreword 

ISO (the International Organization for Standardization) is a worldwide 
federation of national standards bodies (ISO member bodies). The work of 
preparing International Standards is normally carried out through ISO 
technical committees. Each member body interested in a subject for which 
a technical committee has been established has the right to be represented 
on that committee. International organizations, governmental and non- 
governmental, in liaison with ISO, also take part in the work. ISO 
collaborates closely with the International Electrotechnical Commission 
(lEC) on all matters of electrotechnical standardization. 

Draft International Standards adopted by the technical committees are 
circulated to the member bodies for voting. Publication as an International 
Standard requires approval by at least 75 % of the member bodies casting 
a vote. 

Intemational Standard ISO 12213-1 was prepared by Technical Committee 
ISO/TC 193, Natural gas. Subcommittee SO 1. Analysis of natural gas, 

ISO 12213 consists of the following parts, under the general title Natural 
gas -^Calculation of compression factor 

— Part 1: Introduction and guidelines 

— Part 2: Calculation using molar-composition analysis 

— Pari 3: Catoulation using physical properties 

Annex A forms an integral part of this part of ISO 12213. Annexes B and C 
are for information only. 
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Natural gas — Calculation of compression factor — 
Parti: 

\ Introduction and guidelines 

1 Scope 

This tnternationat Standard specifies methods for the calculation of compression factors of natural gases, natural 
gases containing a synthetic admixture and similar mixtures at conditions under which the mixture can exist only as 
tj agas. 

C " ' 

01 

S The standard is in three parts: part 1 gives an introduction and provides guidelines for the methods of calculation 
I described in parts 2 and 3. 

v 

5 Part 2 gives a method for use where the detailed molar composition of the gas is known. Part 3 gives a method for 

1 use where a less detailed analysis, <:omprising superior calorific value (volumetric basis), relative density, carbon 

2 dioxide content and (if non-zero) hydrogen content, is available. 

o 

f Both methods are applicable to dry gases of pipeline quality within the range of conditions under which transmission 
£ and distribution, including metering for custody transfer or other accounting purposes, are normally carried out. In 
^ general, such operations take place at temperatures between about 263 K and 338 K (approximately - 10 'C to 
I 65 *C) and pressures not exceeding 12MPa (120 bar). Within this range, the uncertainty of prediction of both 
% methods is about ± 0,1 % provided that the input data, including the relevant pressure and temperature, have no 
■§ uncertainty. 

NOTE — Pipeline quality gas is used in this International Standard as a concise term for gas which has been processed so as 
to be suitable for use as industhal, commercial or domestic fuel. Although there Is no formal international agreement upon the 
composition and properties of a gas which complies with this concept, some quantitative guidance is provided in 5.1.1. A 
detailed gas quality specification is usually a matter for contractual arrangements between buyer and seller. 

The method given in part 2 is also applicable (with increased uncertainty) to broader categories of natural gas. 
including wet or sour gases, within a wider range of temperatures and to higher pressures, for example at resen^oir 
or underground storage conditions or for vehicular (NGV) applications. 

The method given in part 3 is applicable to gases with a higher content of nitrogen, carbon dioxide or ethane than 
normally found in pipeline quality gas. The method may also be applied over wider ranges of temperature and 
pressure but with increased uncertainty. 

For the calculation methods described to be valid, the gas must be above its water and hydrocarbon dewpoints at 
the prescribed conditions. 

The standard gives all of the equations and numerical values needed to implement both methods. Verified computer 
programmes are available (see annex B). 



1 
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2 Normative references . 

The following standards contain provisions which, through reference in this text, constitute provisions of this part of 
ISO 12213. At the time of publication, the editions indicated were valid. All standards are subject to revision, and 
parties to agreements based on this part of IS0 12213 are encouraged to investigate the possibility of applying the 
most recent editions of the standards indicated below. Members of lEC and ISO maintain registers of currently valid 
International Standards. 

ISO 6976:1995. Natural gas — Calculation of calorific values, density, relative density and Wobbe index from 
composition. 

ISO 13443:1 996, Natural gas — Standard reference conditions, 

3 Definitions 

For the purposes of the various parts of this International Standard, the following definitions apply. 

3.1 compression factor, Z: The ratio of the volume of an arbitrary mass of gas, at a specified pressure and 
temperature, to the volume of the same mass of gas under the same conditions as calculated from the ideal-gas 
law, as follows: 

2=V^(real)/Vm(ideal) ...(1) 
where Vm(ideal) = RTfp ... (2) 

Thus, Zip, r, y) = pVrr^ip, T, yViRT) ... (3) 

where . 



p is the absolute pressure; 

T is the thermodynamic temperature; 

>• is a set of parameters which uniquely characterizes the gas (in principle, the latter may be the complete 
molar composition or a distinctive set of dependent physico-chemical properties, or a mixture of both); 

Vfjx is the molar volume of the gas; 

R is the molar gas constant in coherent units. 

The compression factor is a dimensionless quantity usually close to unity. 

NOTE — The terms "compressibility factor^ and "Z-faclor' are synonymous with compression factor, 

3-2 density, p: The mass of a given quantity of gas divided by its volume at specified conditions of pressure and 
temperature. 

3.3 molar composition: The term used when the proportion of each component in a homogeneous mixture is 
expressed as a mole (or molar) fraction, or mole (molar) percentage, of the whole. 

Thus the mole fraction ;c,- of component i is the ratio of the number of moles of comp)onenl i in a given volume of a 
mixture to the total number of moles of all the components in the same volume of the mixture. One mole of any 
chemical species is the amount of substance which contains the relative molecular mass in grams. A table of 
recommended values of relative molecular masses is given in ISO 6976. 

For an ideal gas. the mole fraction (or percentage) is identical to the volume fraction (or percentage), but this is not in 
general a sufficiently accurate approximation to real-gas tsehaviour for the purposes of this International Standard. 
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3.4 molar calorific value, H\ The amount of heat which would be released by the complete combustion in air of 
the hydrocarbons in one mole of natural gas in such a way that the pressure at which the reaction takes place 
remains constant and all the products of combustion are returned to the same specified temperature as that of the 
reactants, all of these products being in the gaseous state except for water formed by combustion, which is 
condensed to the liquid state at the specified temperature. 

Note that the molar calorific value only includes the hydrocarbons in the natural gas, i.e. inert components (primarily 
nitrogen, carbon dioxide and helium) and other combustible components (such as hydrogen and carbon monoxide) 
are excluded. 

The specified temperature is 298,1 5 K (25 *C) and the reference pressure is 101 ,325 kPa, 
The term 'molar heating value" is synonymous with "molar calorific value". 

3.5 superior calorific value, Hs (volumetric basis): The amount of heat which would be released by the complete 
combustion in air of all the combustible components in unit volume of natural gas in such a way that the pressure at 
whbh the reaction takes place remains constant and all the products of combustion are returned to the same 
specified temperature as that of the reactants, all of these products being in the gaseous state except for water 
formed by combustion, which is condensed to the liquid state at the specified temperature- 
Note that the superior calorific value includes all the combustible components in the natural gas. . 

The reference temperature at which the volume is measured is 273,15 K (0 ''C) and the specified temperature at 
which combustion takes place is 298.15 K (25 **C). The reference pressure is 101,325 kPa. 

NOTE — Annex D of part 3 of this Intemational Standard gives conversion factors which enable superior calorific values and 
relative densities determined at other reference or specified temperatures, and other reference pressures, including the ISO 
standard reference conditions (see ISO 1 3443). to be used as input data for the calculation method described. 

The terms *gross"» "higher^, "upper^ and "total calorific value" and "heating value" are synonymous with "superior 
calorific value". 

3.6 relative density, di The ratio of the mass of a given volume of natural gas to the mass of dry air of standard 
composition which would be contained in the same volume at the same reference conditions of pressure and 
temperature. 

The relative density includes all the components of the natural gas. 
NOTE — The standard composition of dry air is given in ISO 6976. 

in this International Standard, the reference temperature is 273,15 K (O'C) and the reference pressure is 
101,325 kPa (see note in 3.5), 

The term "specific gravity" is synonymous with "relative density". 

3.7 uncertainty of a predicted compression factor, ± AZ: The range of values Z - AZ to Z + AZ within which the 
(unknown) true value is expected to lie with a confidence level of 95 %. This uncertainty may be expressed either as 
an absolute value or as a percentage. 

Estimates of the 95 % confidence limits are, to the extent that this is practicable, established by comparison of test 
data of low uncertainty with calculated values of Z. 



4 General principles 

The methods recommended use equations which are based on the concept that any natural gas may be uniquely 
characterized for calculation of its volumetric properties either by component analysis or by an appropriate and 
distinctive set of measurable physical properties. These characteristics, together with the pressure and temperature, 
are used as input data for the methods. 
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In the sense that the volumetric behaviour of a gas mixture derives directly from the numbers and types of 
molecular interactions (collisions) which take place, a method which explicitly recognizes each molecular constituent 
of the mixture, and its proportion of the whole, is to some degree more fundamental than alternatives. 

The method given in part 2 of this International Standard uses a detailed molar-composition analysis in which all 
constituents present in amounts exceeding a mole fraction of 0,00005 should be represented. The sum of the mole 
fractions used should be unity to within 0,0001. For a typical distributed (pipeline quality) gas, this includes all 
alkane hydrocarbons up to about Cy or Cg together with nitrogen, carbon dioxide and helium. For gases containing 
a synthetic admixture, hydrogen, carbon monoxide and ethylene are also likely to be significant components. For 
broader categories of gas, other components such as water vapour and hydrogen sulfide need to be taken into 
consideration. 

The equation recommended is known as the AGA8 detailed characterization equation, and will be referred to 
hereafter as the AGA8-92DC equation HI (see annex C). It is a revision of the equation described in AGA 
Report No. 8 12). 

The method given in part 3 of this International Standard uses two distinct physical properties, namely superior 
calorific value and relative density, together with the carijon dioxide content. 

NOTE — In principle, any three from superior calorific value, relative density, carbon dioxide content and nitrogen content 
may be used, the calculation methods being essentially equivalent. However, the set comprising the first three is prefen^ed for 
this International Standard. The reader interested in the use of alternative input variables is referred to the GERGTM5 
documentation (31. 

This method is particularly useful in the common situation where a complete molar composition is not available, but 
may also be preferred for its relative simplicity. For gases containing a synthetic admixture, the amount of hydrogen 
needs to be known. 

The equation recommended is known as the SGERG-88 equation PI. This equation is derived from the MGERG-88 
equation [^1, which uses a detailed molar analysis to characterize the gas. 

The evaluation of both the AGA8-92DC and the SGERG-88 equations has been carried out using a large databank 
of high-accuracy (± 0,1 %) compression factor measurements (most of which are traceable to the relevant 
international metrological standards), compiled for the purpose by the Groupe Europeen de Recherches 
Gazieresfsj and the Gas Research Institute l^l [71. Within the transmission and distribution pressure and 
temperature ranges, the equations are of essentially identical perfomnance. 



5 Guidelines 

5.1 Pipeline quality natural gases 
5.1.1 Pipeline quality gas 

To make a definitive quantitative specification of what does or does not represent pipeline quality natural gas is. for 
several reasons, an elusive and contentious objective. Nevertheless, most transmission and distribution engineers 
are familiar with the general concept and will normally know whether any particular natural gas falls within the scope 
of this term. The following is therefore intended simply to provide summary guidance for any other users, rather than 
to formalize criteria for pipeline quality gas. 

Pipeline quality gas is taken to consist predominantly (mole fraction greater than 0,70) of methane and to have a 
superior calorific value normally within the range 30 MJ-m-3 to 45 MJ-m-s (see table 1). Nitrogen and carbon dioxide 
are the main diluents (each up to a mole fraction of about 0.20). 

Ethane (up to a mole fraction of about 0,10), propane, butanes, pentanes and higher alkanes will usually be present 
in steadily decreasing amounts. Minor amounts of helium, benzene and toluene may be present at mole fractions of 
less than 0.001. For natural gases with a synthetic admixture, hydrogen and carl>on monoxide may be present in 
mole fractions of up to about 0,10 and 0.03, respectively, and there may be small amounts of ethylene. No other 
component, such as those found in wet and sour gases (for example water vapour, hydrogen sulfide or oxygen), is 
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normally present in greater than trace amounts, and there should be no aerosol, liquid or particulate matter present. 
Minor and trace components should be treated as specified in part 2 of this Intemational Standard. 

This way of defining pipeline quality gas is not intended to exclude natural gases of other compositions from being 
transported through pipelines. 

The limits allowable for the purpose of this Intemational Standard are given In table 1 . 

Table 1 — Allowable limits for mole fractions of components 



Component 


Mole fraction 


Main comoonenls 




Methane 


3" 0.70 


Nitrogen 


^ 0,20 


Carbon dioxide 


<S0,20 


Ethane 


^^0,10 


Propane 


0,035 


Butanes 


^0,015 


Pentanes 


^ 0.005 


Hexanes 


^ 0.001 


Heptanes 


^ 0,000 5 


Octanes and above 


« 0.000 5 


Hydrogen 


^0,10 


Carbon monoxide 


-fiO.03 


Helium 


^ 0.005 


Water 


< 0.000 15 


Minor an^ trace <?gn)p<?nent5 




Ethylene ,^ 


0,001 


Benzene 


« 0.000 5 


Toluene 


< 0.000 2 


Argon 


« 0.000 2 


Hydrogen sulfide 


^ 0.000 2 


Oxygen 


^ 0.000 2 


Total unspecified components 


< 0.000 1 



f 5,1.2 Transmission and distribution metering 

^ The main use of this International Standard is expected to be for the calculation of compression factors in 
' applications concerned with transnrussion and distribution of pipeline quality gases. The range of conditions 
encountered in such applications varies from country to country, but almost all will be covered by the limits 

263K^r^338K 

OMPa<p'S 12MPa 

The methods given in parts 2 and 3 apply with equal validity for all conditions within these limits. 



5.1.3 Calculation using a molar-composition analysis 

The AGA8-92DC equation may be used for any pipeline quality gas for which a detailed molar-composition analysis 
is available. The components which the analysis should include are; methane, nitrogen, carbon dioxide, carbon 
monoxide, hydrogen, helium, ethane, propane, butanes, pentanes, hexanes, and {if present at mole fractions 
greater than 0,000 05) higher hydrocarbons up to C-iq. The amount of each minor or trace component specified in 
the lower part of table 1 should be demonstratably within the relevant limit. Any non-negligible amount of a minor or 
trace component should be treated in the manner specified in part 2 of this International Standard. 

Within the ranges quoted in 5.1.2. the calculated compression factor values have the same status (i.e. equal 
validity) as those calculated from superior calorific value, relative density and carbon dioxide content. 
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The method may be used in all applications where the composition is subject to regular or semi-continuous 
determination. 

5.1.4 Calculation using physical properties 

The SGERG-88 equation may be used for any pipeline quality natural gas having a superior calorific value in the 
range 30 MJ*nr3 to 45 MJ nrr3, a relative density in the range 0,55 to 0,80, a known carbon dioxide content and a 
known hydrogen content. 

Within the range quoted in 5.1.2, the calculated compression factor values have the same status (i.e. equal validity) 
as those calculated from a complete molar-composition analysis. 

The method may be used in all applications where and d are subject to regular or continuous detemnination. 

5.1.5 Manufactured gases . 

Neither the AGA8-92DC method given in part 2 nor the SGERG-88 method given in part 3 is specifically intended 
for use with manufactured (synthetic) gases, as these may contain substantial amounts of chemical species which 
are atypical of natural gases, or common species in atypical proportions (see 5.2.3). 

Either method may, however, be used if it can be demonstrated that the composition of the synthetic gas closely 
matches that of a possible true natural gas, all of the components falling within the ranges of concentration given in 
5.1.1. with hydrocarbons above butane (C4) either absent or decreasing regularly with increasing cartran number. In 
the sense that hydrocarbons above C4 are virtually absent, liquefied natural. gas may fall within this category. 

In addition, the SGERG-88 method may be used with natural gases containing a coke-oven gas mixture up to the 
specified concentration limit for hydrogen. The method may not. however, be used with "undiluted" coke-oven gas. 

5.1.6 Predictive uncertainty 

Given exact values of alt of the relevant input variables, the methods given in parts 2 and 3 are both expected to 
achieve an uncertainty of prediction of ±0,1 % in the compression factor for the entire ranges of composition and 
physical properties given in 5.1.1 for pipeline quality gas, and of pressure and temperature quoted in 5.1.2 for 
transmission and distribution applications. 

The only exceptions to this are as follows. For the prediction of Z, by the method given in part 3, for gases 
containing more than a 0,15 mole fraction of nitrogen or 0.05 of carbon dioxide (up to the relevant limits of 0,20), an 
uncertainty within ±0,10% is only maintained up to about lOMPa for nitrogen and 6MPa for carbon dioxide 
(see 5.2.2). 

It is stressed, however, that any uncertainty in the input variables adds further uncertainty to the result. The 
sensitivity of the result to the accuracy of the input variables depends significantly, and in a complicated manner, 
upon: 

a) the magnitude of each input variable; 

b) the degree of independeru:e of each input variable from the values of other input variables. 

In most cases, the greatest sensitivity of the result towards all of the input variables is found at the upper extremity 
of the pressure range (12 MPa) and the lower extremity of the temperature range (263 K). 

As a general guideline only, uncorrelated uncertainties of the variables listed in table 2 may contribute to an 
additional uncertainty of the result of about ± 0,1 % at 6 MPa and within the temperature range 263 K to 338 K. 
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Table 2 — Allowable uncertainties of Input variables for A2< 0,1 % 



Input variable 



Allowable uncertainty 



Calorific value 
Relative density 
Temperature 
Pressure 
Mote fraction of 

inerts 

methane 

ethane 

propane 

butane 

pentanes plus higher hydrocarbons 
hydrogen and carbon monoxide 



±0,06 MJ.m-3 
±0.001 3 
±0.15 K 
± 0.02 MPa 



±0.001 

±0,001 

±0.001 

±0.000 5 

±0.000 3 

±0,0001 

±0.001 



Thus, the choice of calculation method should take into account not only the availability of an appropriate form of 
input data, but also its accuracy. In circumstances where it Is unclear whether the instrumental accuracy is 
sufficient, the user should carry out sample calculations of the compression factor at the highest pressure and 
lowest temperature of interest in order to establish, for gases typical of the application in question, the sensitivity to 
independent small variations of all input variables. 

5.1.7 Wider ranges of pressure and temperature 

Both the AGA8-92DC and SERG-8d equations degrade in accuracy outside the primary ranges of pressure and 
temperature given in 5.1 .2. 

However, the AGA8-92DC equation is generally expected to extrapolate beyond these ranges more accurately than 
the SGERG-88 equation, and consequently would often be the preferred alternative for applications outside the 
normal range of transmission and distribution conditions. 

The predictive uncertainty of the AGAB-92DC equation depends strongly upon both the composition of the gas and 
the conditions under consideration. Further advice is given in part 2 concerning the uncertainty of prediction to be 
expected at any selected conditions. This uncertainty may sometimes be almost as low as that for the transmission 
and distribution region, but at extremes of temperature or pressure will be significantly greater. Because of the lack 
of high-accuracy test data, it is extremely difficult to assess the uncertainty. 



5.2 Other gases and other applications 
5.2.1 Introduction 

Not all fluids with which the gas engineer may typically have to deal, and for which he may need to know 
compression factors, are natural gases of pipeline quality. For exannple. unprocessed (welt-head) or partially 
processed natural gases do not normally fiall within the scope of the term pipeline quality gas as defined in 5.1,1, for 
which this International Standard is primarily intended. Nor do manufactured gases. 

Nevertheless, the methods recommended in parts 2 and 3 can be applied, with certain restrictions and with 
increased uncertainty, to such gases. Although such gases may not usually be distributed to end-users, many of the 
applications for which calculations are required are for pressure and temperature values which fall within the ranges 
specified in 5.1.2. 



Self-evidently. the guidelines which can be given (as well as the calculation uncertainties) become less definitive as 
limitations on the allowable pressure, temperature and composition ranges are relaxed. 
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5St.J2 Lean and rich gases 

Some natural gases exist, and are distributed, which contain nitrogen, carbon dioxide, ethane or higher 
hydrocarbons in amounts which exceed the limits for which an uncertainty of ± 0,1 .% can properly be claimed. In 
this context, gases which contain more than a 0,15 mole fraction of nitrogen or a 0,05 mole fraction of carbon 
dioxide are termed "lean", and gases which contain more than a 0,10 mole fraction of ethane or 0.035 of propane, 
and so on. are termed '*rich\ 

The methods recommended in parts 2 and 3 can both the applied to these types of lean and rich gases, but with 
some increase in uncertainty of prediction. 

For example, the method given in part 3 may be applied to give an uncertainty within ± 0,2 % for natural gases 
containing up to about a 0,50 mole fraction of nitrogen, 0,18 of carbon dioxide or 0,13 of ethane at pressures up to 
10 MPa. A more detailed estimate of the uncertainty for each method as a function of component mole fraction is 
given in parts 2 and 3 where plots are given which show, for a wide range of temperatures, the pressurie versus 
mole fraction surfaces for nitrogen, carbon dioxide, ethane and propane, respectively, with uncertainty of prediction 
as a parameter. The major problem in providing such plots is the paucity of high-accuracy test data. 

5.2.3 Wet and sour gases 

This category of gases is taken to comprise those gases which fall short of qualifying as pipeline quality natural 
gases only by the inclusion of undesirable components. Typically such gases may be unprocessed (well-head) or 
partially processed natural gases and may contain, for example, water vapour ("wef gases), hydrogen sulfide 
("sour' gases) or oxygen in amounts significantly greater than those quoted in 5.1,1, perhaps also with traces of 
carbonyl sulfide and process-fluid vapours such as methanol or glycols. 

The method given in part 2 is applicable to any such gas, provided that the unwanted components are limited to 
water vapour, hydrogen sulfide and oxygen. The uncertainty of prediction is, however, substantially increased. The 
method given in part 3 should not be applied to such gases. 

5.2.4 Manufactured gases 

Manufactured gases fall into two distinct categories, as follows: 

a) those which are intended as synthetic or substitute natural gases, arid which closely match true natural gases 
in both composition and properties; 

b) those which, whether or not intended to replace or enhance natural gas in service, do not closely match natural 
gases In composition. 

In case (a) it is clear that, if the composition is such that the gas is indistinguishable from that of a possible true 
natural gas, then the methods given in parts 2 and 3 apply with no increase in uncertainty (see 5.1.5). This is, 
however, rarely likely to be the case. More often the manufactured gas, even if it contains inert and tower 
hydrocarbons in satisfactory proportions, will not exhibit the distinctive hydrocarbon "tail" of a true natural gas and 
may additionally contain small but significant amounts of non-aikane hydrocarbons. It is difficult to assess the 
effects of this upon the uncertainty of prediction. 

Case (b) includes gases such as town gas, (undiluted) coke-oven gas, and LPG-air mixtures, none of which is 
compositionally similar to a true natural gas (even though, in the latter case, it may be operationally interchangeable 
with natural gas). The method given in part 3 should not be applied to any such gases. The method given in part 2 
may be applied, but the uncertainty of prediction is extremely difficult to assess. 

5.2.5 Summary of predictive uncertainty 

The expected uncertainties of the calculation methods given in part 2 and part 3 are summarized in figure 1. The 
uncertainties are given as pressure versus mole fraction bar charts for various contents of nitrogen, carbon dioxide 
and ethane. Results are summarized for. 

a) pipeline quality natural gases within the normal ranges of temperature and pressure for transmission and 
distribution applications (263 K to 338 K; 0 MPa to 1 2 MPa); 
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b) pipeline quality natural gases (or the same range of temperature, but for a wider range of pressure (up to 
30MPa); 

c) wider ranges of gas composition (up to 0,5 N2, 0,3 CO2 and 0,2 C2H6) for the same ranges of temperature and 
pressure. 

The more detailed inforrhation on which these bar charts are based is given in annex E of parjs 2 and 3. Information 
concerning the performance outside the primary range of temperature and pressure is given in figure 1 of parts 2 
and 3. 




^K7 ^COZ *C02 'fCO? 'C2HA ^C2H4 ^ N2 'CM 'C2HA 

0,0s 0.10 0.20 O.OS 0.10 0.20 0.05 0.10 0.50 0.30 0.20 



^ D AGA8-DC92 calculation method (part 2) using motar-composition analysis 

S SGERG-88 calculation method (part 3) using physical properties 

\) Pipeline quality natural gases (263 K to 338 K; 0 MPa to 12 MPa) 

(D) Pipeline quality natural gases (263 K to 338 K; 12 MPa to 30 MPa) 

(c) Wider ranges of gas composition (263 K to 338 K; 0 MPa to 30 MPa) 

Figure 1 — Expected uncertainty of the calculation methods given in parts 2 and 3 
5.2.6 Calculation of related properties 

Although the express purpose of this International Standard is to permit calculations of compression factors, it is 
appropriate to note that other properties of natural-gas-type fluids may also be calculated by means of the methods 
described in parts 2 and 3. Self-evidently, the molar density p^, being simply the reciprocal of the molar volume 
Vm(real), is always available from equations (1) and (2) once Z{p, T) is known. The (mass) density p is also 
available, as the product of pm and the mean molar mass M (molecular weight), if the latter is known, as is the case 
for any fluid characterized by a molar-composition analysis. 

If the molar composition is not known, the mass density can instead be calculated by using the compression factors 
at pipeline and normal conditions together with the relative density and known mass density of dry air at normal 
conditions (see part 3, equation (8.42)). 
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Annex A 

(normative) 

Symbols and units 



Symbol 


Meaning 


d 


Relative density 


H 


Molar calorific value 


Hs 


Superior calorific value 


M 


Molar mass 


p 


Absolute pressure 


R 


Gas constant 


T 


Absolute temperature 


Vm 


Molar volume 




Molar (mole) fraction of component / 


y 


. Set of properties 


z 


Compression factor 


bZ 


Uncertainty (95 % confidence fimits) of predicted 




compression factor 


9 


Mass density 


Pm 


Molar density 



Value 


Units 


VanaOie 




variable 


kJ-mol-i 


variable 


MJ.m-3 


variable 


kg-kmol-i 


variable 


kPa 


8.314 510 


JmoI-i.K-i 


variable 


K 


variable 


m3.kmol-i 


variable 




variable 




variable 




variable 


kg-nrr3 




kmolnrr3 
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Annex B 

(informative) 

Suppliers of computer programmes 

it is planned to malce software available wliich implements tliis International Standard. Users are invited to contact 
their ISO member body or ISO Cientral Secretariat to enquire about the availability of such software. 
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(informative) 
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Foreword 

ISO (the International Organization for Standardization) is a worldwide 
federation of national standards bodies (ISO nnember bodies). The work of 
preparing International Standards is nornnally carried out through ISO 
technical committees. Each member body interested in a subject for which 
a technical committee has been established has the right to be represented 
on that committee. International organizations, governmental and non- 
governmental, in liaison with ISO, also take part in the work. ISO 
collaborates closely with the International Electrotechnical Commission 
(lEC) on all matters of electrotechnical standardization. 

Draft International Standards adopted by the technical committees are 
circulated to the member bodies for voting. Publication as an International 
Standard requires approval by at least 75 % of the member bodies casting 
a vote. 

International Standard ISO 12213-2 was prepared by Technical Committee 
ISO/TC 193, Natural gas, Subcommittee SC 1, Analysis of natural gas, 

ISO 12213 consists of the following parts, under the general title Natural 
gas — Calculation of compression factor 

— Parf 1: Introduction and guidelines 

— Part 2: Calculation using molar-composition analysis 

— Part 3: Calculation using physical properties 

Annexes A to D form an integral part of this part of ISO 12213. Annexes E 
to G are for information only. 
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Natural gas — Calculation of compression factor — 
Part 2: 

Calculation using molar-composition analysis 

1 Scope 

This International Standard specifies methods for the calculation of compression factors of natural gases, natural 
gases containing a synthetic admixture and similar mixtures at conditions under which the mixture can exist only as 
a gas. 

This part of ISO 12213 specifies a method for the calculation of compression factors when the detailed composition 
of the gas by mole fractions is known, together with the relevant pressures and temperatures. 

The method is applicable to pipeline quality gases within the ranges of pressure p and temperature T at which 
transmission and distribution operations normally take place, with an uncertainty of about ± 0.1 %. It can be applied, 
with greater uncertainty, to wider ranges of gas composition, pressure and temperature (see annex E). 

More detail concerning the scope and field of application of the method is given in part 1 of this International 
Standard. 



2 Normative references 

The following standards contain provisions which, through reference in this text, constitute provisions of this part of 
ISO 12213. At the time of publication, the editions indicated were valid. All standards are subject to revision, and 
parties to agreements based on this part of ISO 12213 are encouraged to investigate the possibility of applying the 
most recent editions of the standards indicated below. Members of lEC and ISO maintain registers of currently valid 
International Standards. 

ISO 31 -3:1 992, Quantities and units — Part 3: Mechanics, 
ISO 31-4:1992, Quantities and units — Part 4: Heat, 

150 6976:1995, Natural gas — Calculation of calorific values, density, relative density and Wobbe index from 
composition, 

ISO 12213-1:1997, Natural gas — Calculation of compression factor — Part 1: Introduction and guidelines. 



3 Definitions 

All definitions relevant to the use of this part of ISO 12213 are given in part 1 . 
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4 Method of calculation 

4.1 Principle 

The method recommended uses an equation based on the concept that pipeline quality natural gas may be 
uniquely characterized for calculation of its volumetric properties by component analysis. This analysis, together 
with the pressure and temperature, are used as input data for the method. 

The method uses a detailed molar-composition analysis in which ail constituents present in amounts exceeding a 
mole fraction of 0,000 05 should be represented. Typically, this includes all alkane hydrocarbons up to atjout C7 or 
Ca together with nitrogen, carbon dioxide and helium 

For other gases, additional components such as water vapour, hydrogen sulfide and ethylene need to be taken into 
consideration (see reference [1] in annex G). 

For manufactured gases, hydrogen and carbon monoxide are also likely to be significant components, 

4.2 The AG A8-92DC equation 

The compression factor is determined using the AGA8 detailed characterization equation (denoted hereafter as the 
AGA8-92DC equation). This is an extended virial-type equation. The equation is described in AGA Report No. Sl^l, It 
may be written as 

Z = 1+ SAn-Pr £c;+ 2c(^n-cA/V^«)p^exp(-C,p^ ' . . . (1) 

n = 13 /I a 13 



Is the compression factor; 
is the second virial coefficient; 
is the molar density (moles per unit volume): 
is the reduced density; 
, are constants (see table B.1); 
are coefficients which are functions of temperature and composition. 

The reduced density pr is related to the molar density pm by the equation 

Pr'K^Pm •••(2) 
where ^is a mixture size parameter. 
The molar density can be written as 

p^^p/{ZJtT) ...(3) 

where 

p is the absolute pressure; 

R is the universal gas constant; 

T is the absolute temperature. 

Z is calculated as follows: first the values of B and C* (n = 13 to 58) are calculate^, using relationships given in 
annex B. Equations (1) and (3) are then solved simultaneously for pm and Z by a suitable numerical method (see 
figure 8.1). 
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Z = 


c 




1 


where 
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Z 
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U 
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Pm 
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C* 
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4.3 Input variables 

The input variables required for use witli the AGA8-92DC equation are the absolute pressure, absolute temperature 
and molar composition. 

The composition is required, by mole fraction, of the following components: nitrogen, carbon dioxide, argon, 
methane, ethane, propane, n-butane, methyl-2-propane (iso-butane), /t-pentane, methyl-2-butane (iso-pentane), 
hexanes. heptanes, octanes, nonanes. decanes, hydrogen, carbon monoxide, hydrogen sulfide, helium, oxygen 
and water. 

NOTE — If the mole fractions of the heptanes, octanes, nonanes and decanes are unknown, then use of a composite C^^ 
fraction may be acceptable. The user should carry out a sensitivity analysts in order to test whether a particular approximation 
of this type degrades the result. ' 

All components with mole fractions greater than 0,000 05 shall be accounted for. Trace components (such as 
ethylene) shall be treated as given in table 1. . 

If the composition is known by volume fractions, these shall be converted to nniole fractions using the method given 
in ISO 6976. The sum of all mole fractions shall be unity to within 0,000 1. 

4.4 Ranges of application 
4.4.1 Pipeline quality gas 

The ranges of application for pipeline quality gas are as defined below: 

absolute pressure 0 MPa ^12MPa 

temperature 263 K ^ T ^ 338 K 

superior calorific value 30 MJ-nrt^ ^ Ws ^ 45 MJ-m^ 

relative density 0,55 ^ d ^ 0,80 

The mole fractions of the natural-gas components shall be within the following ranges: 



methane 


0,7 


' ^ J^CH4 


^ 1,00 


nitrogen 


0 




<0,20 


carbon dioxide 


0 




=5 0,20 


ethane 


0 




^0,10' 


propane. 


0 


J^caHa 


^ 0.035 


butanes 


0 


'C4H10 


=s 0,015 


pentanes 


0 




< 0,005 


hexanes 


0 




^ 0,001 


heptanes 


0 




^ 0.000 5 


octanes plus 


0 




^ 0.000 5 



higher hydrocarbons 

hydrogen 0 jch2 '^0,10 

carbon monoxide 0 « xqo ^ 

helium 0 XHe ^ 0,005 

water 0 ^ XH2O ^ 0.000 1 5 

Any component for which x,- is less than 0.000 05 can be neglected. 
Minor and trace components are listed in table 1. 
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Table 1 — Minor and trace components 







Oxygen 


oxygen 


Argon 


argon 


Hydrogen sulfide 


hydrogen sulfide 


Ethylene, acetylene 


cart>on dioxide 


Propylene, propadiene 


propane 


134*4 ARAM t#^^AM*% 

DUienes, Duiaaienes 


n-butane 


Neo-pentane, pentenes, benzene, cyctopentane . 


n-pentane 


Alt Cs-Isomers, cyctohexane, ethylbenzene, xylenes 


/t-hexane 


All C7-fsomers, cycloheptane, toluene 


n-heptane 


AnCg-isomers 


. n-octane 


All Cg-isomers 


n-nonane 


All C^o'l^'^®''^ higher hydrocarbons 


n-decane 



The method applies only to mixtures in the single-phase gaseous state (above the dew point) at the conditions of 
temperature and pressure of interest. 

4.4.2 Wider ranges of application 

The ranges of application tested beyond the limits given in 4.4.1 are: 



absolute pressure 
temperature 
relative density 
superior caloriftc value 



0 MPa ^65 MPa 

225 K «r ^350K 
0,55 . ^ 0.90 

20 MJ-mri ^ Hs « 48 MJ nrS 



The allowable mole fractions of the major natural gas components are: 



methane 

nitrogen 

carbon dioxide 

ethane 

propane 

hydrogen 



0,50 'SJCCH4 ^ 1»00 
0 jcr42 ^ ^'5^ 
0 *5jtco2 ^0.30 
0 =^JccaH6'S^20 

0 ^jcHz ^0,10 



The limits for minor and trace gas components are as given in 4.4.1 for pipeline quality gas. For use of the method 
outside these ranges, see annex E. 



4.5 Uncertainty 

4.5.1 Uncertainty for pipeline quality gas 

The uncertainty of results for use on all pipeline quality gas within the limits described in 4.4.1 is ±0.1 % (for the 
temperature range 263 K to 350 K and pressures up to 12 MPa) (see figure 1). For temperatures above 290 K and 
at pressures up to 30 MPa the uncertainty of the result is also ± 0,1 %. 



For lower temperatures, the uncertainty of ± 0,1 % is at least maintained for pressures up to about 10 MPa. 
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AGA6^92 equation 



A a sO.1 % 
' B £0.1% to ±0,2% 
C *0.2 % to £0,S % 



250 275 300 325 350 

Tempepature. r(K) 

Figure 1 — Uncertainty limits for the calculation of compression factors (The uncertainty limits given are 
expected to be valid for natural gases and similar gases with ^ ^COz 0-20, .VC2H6 ^ 0.^0 
and JCH2 < 0,10, and for 30 MJm-3 ^Hs'^^S MJm-3 and 0,55 0,80) 



This uncertainty level has been determined by comparison with the GERG databank of measurements of the 
compression factor for natural gases I^]. 13]. A detailed comparison was also made with the GRI PVT data on 
gravimetrically prepared simulated natural-gas mixtures W» [^1. 

The uncertainty of the measurements in both databanks used to test the method is of the order of ± 0,1 %. 
4^.2 Uncertainty for wider ranges of application 

The estimated uncertainties for calculations of compression factors beyond. the limits of quality given in 4.4.1 are 
discussed in annex E. 



4,5.3 Impact of uncertainties of input variables 

Listed in table 2 are typical values for the uncertainties of the relevant input variables. These values may be 
achieved under optimum operating conditions. 

As a general guideline only, an error propagation analysis using the uncertainties in the input variables produces an 
additional uncertainty of about ± 0.1 % in the result at 6 MPa and within the temperature range 263 K to 338 K. 
Above 6 MPa. the additional uncertainties are greater and increase roughly in direct proportion to the pressure. 
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Table 2 — Uncertainties of input variables 



Input variable 


Absolute uncertainty 


Absolute pressure 


±0,02MPa 


Temperature 


±0J5K 


Mole fraction of 




inerts 


•4-onni 


nitrogen 


±0,001 


ceufbon dioxide 


±0.001 


methane 


±0,001 


ethane 


±0.001 


propane 


±0.000 5 


butanes 


±0.000 3 


pentanes plus higher hydrocarbons 


± 0,000 1 


hydrogen and carbon monoxide 


±0.001 



4,5.4 Reporting of results. 

Results for compression factor and molar density shall be. reported to four and to five places of decimals, 
respectively, together with the pressure and temperature values and the calculation method used (ISO 12213-2, 
AGA8-92DC equation). For verification of calculation procedures, it is useful to carry extra digits, 

5 Suppliers of computer programmes 

It is planned to make software available which implements this International Standard, Users are invited to contact 
their ISO member body or ISO Central Secretariat to enquire about the availability of such software. 
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Annex A 

(normative) 

Symbols and units 



c 

3 



e 



Symbol 

On 
B 

B„ 

E) 

F 

Fi 

Pj 
f» 
C 

Ci 

Gj 

gn 

Hs 

K 



M 

Mi 

N 

n 

P 

Q 
Qi 
Qj 

R 

Si 



Meaning 

Constant in table B.I 
Second virial coefficient 

Mixture interaction coefficient [equations (B.1) and (B.2)] 

Constant in table B.I 
Constant in table B.I 

Coefficients which are functions of temperature and composition 

Characteristic energy parameter for ith component (table B.2) 

Characteristic energy parameter for;th component (table B.2) 

Binary energy parameter for second virial coefficient 

Binary energy interaction parameter for second virial coefficient (table B.3) 

Mixture high-temperature parameter 

High-temperature parameter for ith component (table B.2) 

High-temperature parameter for j\h component (table B.2) 

Constant in table B,1 

Mixture orientation parameter 

Orientation parameter for ith component (table B.2) 

Orientation parameter for Jth component (table B.2). 

Binary orientation parameter 

Binary interaction parameter for orientation (table B.2) 

Constant in table B.I 
Superior calorific value 
Size parameter 

Size parameter for ith component (table B.2) 
Size parameter for ^th component (table B.2) 
Binary interaction parameter for size (table B.3) 
Constant in table B.I 
Molar mass 

Molar mass of ith component 

Number of components in gas mixture 

An integer (from 1 to 58) 

Absolute pressure 

Quadrupole parameter 

Quadrupole parameter for ith component 

Quadrupole parameter for j\h component 

Constant (table B.I) 

Gas constant (= 0,008 314 510) 

Dipole parameter for ith component (table B.2) 

Dipole parameter for jih component (table B.2) 

Constant (table B.I) 



Units 

m3kmol-i 



K 
K 
K 



MJ.m3 

(m3/kmol)i/3 
(m3/kmol)^/3 

(m3/kmol)i/3 



kgkmol-1 
kg-kmol-1 



MPa 



MJ.(kmolK)- 
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Symbol Meaning 

T Absolute temperature 

V Mixture energy parameter 

Uij Binary interaction parameter for mixture energy (table B.3) 

u„ Constant in table B.I 

Wi Association parameter for ilh component (table B.2) 

• Wj Association parameter for yth component (table B.2) 

Wn Constant (table B.1) 

Xi Mote fraction of ilh component in gas mixture 

xj Mole fraction of 7th component in gas mixture 

Z Compression factor 

p Mass density 

Pr Reduced density of gas 

Pm Molar density 
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Annex B 

(normative) 



Description of the AGA8-92DC method 



B.I General 

For gas mixtures, the compression factor Z is calculated using the equations given in 4.2. This annex gives a 
detailed description of the computations and the necessary numerical values. The description is based upon that 
given in AGA Report No. 8 PI. A programme implementing this descriptiori is given in annex F, and as such 
provides the correct solution. Other computational procedures are acceptable provided that they can be 
demonstrated to yield identical numerical results (see annex O for examples). 



B.2 Computer implementation of the AGA8-92DC method 
B.2.1 Overview of the calculation procedure 

1 Input the absolute temperature 7, absolute pressure p and mole fraction of each component of the mixture. 

NOTE — For pressure and temperature, values known in any other units will first have to be converted precisely to 
values in megapascals and Kelvins, respectively (see ISO 31-3 and ISO 31*4 and annex D for relevant conversion 



II Compute the equation of state coefficients B and Cji (n = 13 to 58) that depend on T and x^. 

III Solve iteratively for the molar density pm. using the equation of state rearranged to give the pressure p, 

IV Output the compression factor after the computed pressure from step III and the input pressure from step 1 
agree within a specified range of convergence (e.g. 1 E-06). 

Figure B.I shows a flow diagram of these steps. 
B.2.2 Details of the calculation procedure 



Input the absolute temperature T, absolute pressure p and mole fraction Xi of each constituent in the natural-gas 
mixture. 



At the absolute temperature r and the mole fractions x,- of the natural gas (as input from step!), compute the 
composition- and temperature-dependent coefficients B and (n = 13 to 58). 

The second viria! coefficient B is given by the following equations: 



factors). 



Step I 



Step II 




...(B.I) 



. . . (B.2) 
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Input T, p. 



Compute 
equation of state coefficients 
BzntiCn tn- 13 to 58) 



Solve 
(density searcti) 




Yes 



Output 



No 



Rgure B.I — AGA8-920C equation — Calculation flow diagram 



The binary parameters E-^ and G,y are calculated using the following equations: 

£i;=£^(£i£^)^ ...(B.3) 

Cii = C'ij(Ci + C;)/2 ... (B.4) 

Note that all values of the binary interaction parameters £,y and are 1,0 except for the values given in table B.3. 
The coefficients C'„ (n = 13 to 58) are given by the equation: 

c;=a„(C + 1-«„)*'(62+1-9„)'-(F + 1-/,/"t;''«r-«« ...(B.S) 

The mixture parameters U, C and Q are calculated using the following conformal solution mixing equations, where 
in the double sums i ranges from 1 to N - 1 and, for each value of i, j ranges from i + MoN: 



ty5 = 



0 = 1 



1 = 1 



1=1 y a I + 1 
/asl y = i + i 



. . . (B.6) 



(B.7) 
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. . . (B.8) 



i>1 



...(B.9) 



It should be noted that all values of the binary interaction parameters Kij, Ejj, Cjj and Uij are 1,0 except for the 
values given in table B.3. Also note that Fi is zero for all components except hydrogen, for which FiHz) - 1,0. and 
that Wi is zero for all components except water, for which W(H20) = 1 ,0, 

Step III 

In the computation of the compression factor Z, the composition of the gas is known, the absolute temperature T of 
the gas is known and the absolute pressure is known. The problem then is to compute the molar density pr^, using 
the equation of state expression for the pressure p. For this purpose, the definition of the compression factor Z as 
given in equation (1) (see 42) is substituted into equation (3) to obtain an equation for the pressure as given in 
equation (B.10): 



18 58 



n = 13 



n«=13 



• . . . (8.10) 



Equation (8.10) is solved using standard equation of state density search algorithms. Having obtained an equation 
for the pressure p [equation (8.10)]. the problem is then to search for the value of the molar density pm that will yield 
the pressure that is within a preset limit (e.g. 1 x ^(H) equal to the input pressure. 

The reduced density pf is related to the molar density pm by the mixture size parameter [see equation (2) in 4.2]. 
The mixture size parameter K is calculated using the following equation: 



a:5 = 



/V-1 



N 



+ 2£ .Xx,x,.(4-i)(/r,/f,) 

,=1 



...(B.11) 



Note that in the summations the subscript i refers to the ith component in the gas mixture and the subscript J refers 
to the yth component in the gas mixture. The quantity N is the number of components in the mixture. Thus, in the 
single summation, / ranges over the integer values from 1 to N, For example, for a mixture of 12 components, 
N = 12 and there would be 12 terms in the single sum. in the double summation, i ranges from 1 to N-^ and, for 
each value of /,/ ranges from i + 1 to M For example, for a mixture of 12 components, there would be 66 terms in 
the double summation if all values of Kij differed from 1,0. However, because n^ny of the values of Ky are 1,0, the 
number of non-zero terms in the double sumnr^tion is small for many natural-gas mixtures. Note that all values of 
Kij are 1 ,0 except for the values given in table 8.3. 



Step IV 

Once the molar density pm has been obtained in step III, the compression factor is calculated in step IV using the 
pressure, temperature, molar density and gas constant: 



Z^p/if^RT) 



NOTE — The density p (mass per unit volumie) can be calculated as follows: 



.(B.12) 
.(8.13) 
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where Af is calculated from the equation: 
1 = 1. 

Report the density to three places of decimals. 



1 0,153 832 600 

2 1,341 953 000 

3 - 2.998 583 000 

4 -0.048 312 280 

5 0,375 796 500 . 

6 -1,589 575 000 

7 -0.053 588 470 

6 0.886 594 630 . 

9 -0.710 237 040 

10 -1.471 722 000 

11 1.321 650 350 

12 -0.786 659 250 

13 2,291 290x10-9 

14 0.157 672400 

15 -0.436 386 400 

16 -0,044 081 590 

17 -0,003 433 888 

18 0,032 059 050 

19 0.024 873 550 

20 0.073 322 790 

21 - 0.001 600 573 

22 0.642 470 600 

23 -0.416 260 100 

24 -0.066 899 570 

25 0.279179 500 

26 -0.696 605 100 

27 - 0,002 860 589 

28 - 0.008 098 836 

29 3.150 547 000 

30 0,007 224 479 

31 -0.705 752 900 

32 0.534 979 200 

33 - 0.079 314 910 

34 -1.418 465 000 

35 -5.999 05x10-^7 

36 0.105 840 200 

37 0.034 317 290 

38 - 0.007 022 847 

39 0.024 955 870 

40 0.042 968 180 

41 0,746 545 300 

42 - 0,291 961 300 

43 7.294 616 000 

44 -9.936 757 000 

45 -0.005 399 808 

46 -0.243 256 700 

47 0.049 870 160 

48 0.003 733 797 

49 1,874 951 000 

50 0.002 168 144 

51 -0,658 716 400 

52 0,000 205 518 

53 0.009 776 1 95 

54 - 0.020 487 080 

55 0,015 573 220 

56 0.006 862 415 

57 -0.001 226 752 

58 0.002 850 908 



Table B.1 — Equation of state parameters 



0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

2 

2 

2 

4 

4 

0 

0 

2 

2 

2 

4 

4 

4 

4 

0 

1 

1 

2 

2 

3 

3 

4 

4 

4 

0 

0 

2 

2 

2 

4 

4 

0 

2 

2 

4 

4 

0 

2 

0 

2 

1 

2 

2 

2 

2 



0.0 
0.5 
1.0 
3,5 
-0.5 
4,5 
0.5 
7,5 
9,5 
6,0 
12,0 
12.5 
-6.0 
2.0 
3.0 
2.0 
2.0 
11.0 
-0.5 
0.5 
0,0 
4.0 
6.0 
21 .0 
23,0 
22.0 
-1.0 
-0.5 
7,0 
-1,0 
6.0 
4.0 
1.0 
9.0 
-13.0 
21,0 
8.0 
-0.5 
0.0 
2.0 
7.0 
9.0 
22,0 
23.0 
1,0 
9.0 
3.0 
8.0 
23,0 
1.5 
5.0 
-0.5 
4.0 
7.0 
3.0 
0.0 
1,0 
0.0 



0 
0 
0 
0 

1 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 
0 

1 

0 . 
0 

1 
1 
1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

6 

0 

1 

0 
0 

1 

0 

1 

0 
0 



^« 

0 
0 
0 
0 
0 
0 

1 

0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 

1 

0 

1 

0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 
0 
0 

1 

0 
0 

. 0 
0 

1 

0 

1 

0 
0 

1 

0 
0 
0 
0 
0 
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Table B*2 — Characterization parameters 



IdenUfi- 


Compound 


Malar mafitt 




9i<C9 par dineujr 


vnvnuiuon 


uuaorupoiB 


HI ah* 


Dipole 


Assocl- 


cation 












nam mAf er 
pal aincivv 




p<ii Miiivwr 




number 












param- 




parafTt- 
















eter 




eter 






^fi 








Ot 














K 


{m3/kmoI)^'3 












1 


Methane 


16.043 0 


151,318 300 


0.461 925 5 


0,0 


0.0 


0,0 


0,0 


0,0 


2 


Nitrogen 


' 28.013 5 


99,737 780 


0,447 915 3 


0.027 815 


0,0 


.0.0 


0,0 


0.0 


3 


Carbon dioxide 


44,010 0 


241.960 600 


0,455 746 9 


0.169 065 


0.690 000 


0,0 


0,0 


0,0 


4 


Ethane 


30,070 0 


244,166 700 


0,527 920 9 


0,079 300 


0.0 


0,0 


0,0 


0 0 


5 


Propane 


44.097 0 


298.118 300 


0,583 749 0 


0,141 239 


0.0 


0.0 


0,6 


00 


6 


Water 


18 0153 


514 015 600 


0,382 586 8 * 


0 332 500 


1.067 750 


0.0 


1 ,582 200 


1 0 


7 


Hydrogen sulfide 


34,082 0 


296 355 000 


0,461 826 3 


0,088 500 


0,633 276 


0,0 


0.390 000 

W w w w w\/ 


0,0 


8 


Hydrogen - 


2,015 9 


26 957 940 


0,351 491 6 


0034 369 

WpWW^T W WW 


00 


1.0 


0.0 


0.0 


9 


Carbon monoxide 


28 0100 


1 05 534 800 


0 453 389 4 


0 038 953 

Wf www WWW 


0,0 


0.0 


0.0 


00 


10 


Oxygen 


31 ,998 8 


122.766 700 


0,418 695 4 


0.021 000 


0.0 


0,0 


0,0 


0 0 


11 


iso-Butane 


58.123 0 


324^068 900 


0,640 693 7 


0.256 692 


0.0 


0.0 


0,0 


00 


12 


n-Butane 


58.123 0 


337,638 900 


0,634 142 3 


0,281 835 


0,0 


0,0 


0,0 


0.0 


13 


iso-Penlane 


72.150 0 


365.599 900 


0,673 857 7 


0.332 267 


0.0 


0,0 


0,0 


0.0 


14 


n-Pentane 


72,150 0 


370,682 300 


0.679 830 7 


0,366 911 


0.0 


0,0 


0,0 


0.0 


15 


n-Hexane 


86,177 0 


402,636 293 


0.717 511 8 


0.289 731 


0.0 


0.0 


0,0 


0.0 


16 


n-Heptane 


100,204 0 


427.722 630 


0,752 518 9 


0,337 542 


0.0 


0.0 


0.0 


0,0 


17 


n-Octane 


114,231 0 


450.325 022 


0,784 955 0 


0.383 381 


0.0 


0.0 


0,0 


0.0 


18 


n-Nonane 


128.258 0 


470.840 891 


0.815 273 1 


0.427 354 


0,0 


0,0 


0.0 


0.0 


19 


n-Decane 


142,285 0 


489,558 373 


0,843 782 6 


0.469 659 


0,0 


0,0 


0,0 


0.0 


20 


Helium 


4,002 6 


2.610 111 


0,358 988 8 


0.0 


0.0 


0,0 


0.0 


0,0 


21 


Argon 


39.948 0 


119,629 900 


0.421 655 1 


0.0 


0.0 


0,0 


0.0 


0,0 



Table B.3 — Binary interaction parameter values 



Identification Component pair Ejj Uij Kij O^j 

number 



2 


Methanes nitrogen 


0.971 640 


0,886 106 


1.003 630 


3 


carbon dioxide 


0,960 644 


0,963 827 


0.995 933 


4 


ethane 








5 


propane 


0,994 635 


0,990 877 


1.007 619 


6 


water 


0.708 218 • 






7 


hydrogen sulfide 


0,931 484 


0,736 833 


1.000 080 


8 


hydrogen 


1.170 520 


1,156 390 


1,023 260 


9 


carbon monoxide 


0,990 126 






10 


oxygen 








11 


iso-butane 


1,019 530 






12 


n*butane 


0.989 844 


0.992 291 


0,997 596 


13 


iso-pentane 


1.002 350 






14 


n-pentane 


0.999 268 


1,003 670 


1,002 529 


15 


n-hexane 


1.107 274 


1.302 576 


0,982 962 


16 


n-heptane 


0.880 880 


1,191 904 


0,983 565 


17 


n-octane 


0,880 973 


1.205 769 


0,982 707 


18 


n-nonane 


0,881 067 


1,219 634 


0,981 849 


19 


n-decane 


0,881 161 


1,233 498 


0,980 991 



0.807 653 



1.957 310 
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Identification 
number 



Component pair 



I 

2 



1 



I 

< 



5 
7 



J 

. 3 
4 
'5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 

■ 18 
19 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
8 
12 
15 
16 
17 
18 
19 
9 

10 
11 
12 



Nitrogen + carbon dioxide 
ethane 
propane 
water 

hydrogen sulfide 
hydrogen 
carbon monoxide 
oxygen 
iso-butane 
A*butane 
iso-pentane 
/i-pentane 
Carbon dioxide -i- ethane 
propane 
water 

hydrogen sulfide 
hydrogen 
carbon monoxide 
oxygen 
iso-butane 
n-butane 
iso-pentane 
n-pentane 
n-hexane 
/i-heptane 
n-octane 
n-nonane 
n*decane 
Ethane + propane 
water 

hydrogen sulfide 
hydrogen 
carbon monoxide 
oxygen 
iso-butane 
/i-butane 
iso-pentane 
/i-penlane 
Propane -H hydrogen 
/i-butane 

Hydrogen sulfide + n-hexane 
n-heptane 
n-octane • 
/i-nonane 
n-decane 
Hydrogen + carbon monoxide 
oxygen 
iso-butane 
/t-butane 



V 






1,022 740 


0,835 058 ' 


0,982 361 


0,970 120 


0.816 431 


1,007 960 


0.945 939 


0,915 502 


0,746 954 






0,902 271 


0,993 476 


0,942 596 


1.086 320 


0,408 838 


1,032 270 


1.005 710 




1,021 000 






0.946 914 






0,973 384 






0,959 340 






0.945 520 






0,925 053 


0.969 870 


1.008 510 


0.960 237 




0,849 408 






U,930 052 


1.045 290 


1,007 790 


1,^01 790 




1,500 000 


0.900 000 




0,906 849 






0,897 362 






0.726 255 




■ 


0,859 764 








1,066 638 


0,910 183 


0,831 229 


1.077 634 


0.895 362 


0,808 310 


1.088 178 


0.881 152 


0,786 323 


1.098 291 


0,867 520 


0,765 171 


1.108 021 


0.854 406 


1,022 560 


1.065 173 


0.986 893 


0.693168 




0,946 871 


0,971 926 


0.999 969 


1,164 460 


1.616 660 


1.020 340 



1.250 000 



1.013 060 


1.250 000 






1,250 000 




1.005 320 


1,250 000 




1,034 787 






1.004 900 






1,008 692 


1,028 973 


0.968130 


1,010 126 


1.033 754 


0.962870 


1,011 501 


1,038 338 


0,957828 


1,012 821 


1,042 735 


0.952 441 


1.014 089 


1.046 966 


0.948 338 


1.100 000 




1.300 000 






1,300 000 







0.982 746 



0,370 296 
1,673 090 
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Annexe 

(normative) 

Example calculations 

The following example calculations have been carried out using the validated computer programme described in 
reference [1], which incorporates the subroutine described in annex B. 



Table C.I — Gas analysis in molis fractions 





Gas1 


Gas 2 


Gas 3 


Gas 4 


Gas5 


Gas 6* 






0.006 


0,005 


0.015 • 


0,016 


0,076 


0,011 






0.003 


0.031 


0.010 


0,100 


0,057 


0.117 


A*H2 




0,00 . 


0.00 


0.00 


0,095 


0.00 


0.00 


XCO 




0.00 


0.00 


0,00 


0,010 . 


0.00. 


0,00 






0,965 


0,907 


0,859 


0,735 


0.812 


0,826 






0.018 


0,045 0 


0.085 


0,033 


0,043 


0.035 






0.004 5 


0.008 4 


0,023 


0.007 4 


0,009 


0,007 5 


''^iso-C4Hio- 


0.001 0 


0.001 0 


0.003 5 


0.001 2 


0,001 5 


0.0012 






0.001 0 


0.001 5 


0,003 5 


0,001 2 


0,001 5 


0.001 2 


•'^iso-CgHia 


0.000 5 


0,000 3 


0.000 5 


0.000 4 


0.00 • 


0.000 4 






0,000 3 


0,000 4 


0,000 5 


0.000 4 


0,00 


0.000 4 


JfC6Hi4 




0.000 7 


0.000 4 


0,00 


0.000 2 


0,00 


0.000 2 






0,00 


0,00 


0,00 


0,000 1 


0.00 


0.000 1 






0.00 


0.00 


0.00 


0,000 1 


0.00 


0,00 


Table — Results (Z-values) 


Cond 


lUons 


Ga$ 1 


Gas 2 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


P 

bar 


•c 














60 


-3.15 


0,840 53 


0.833 48 


0,793 80 


0.885 50 


0.826 09 


0,853 80 


60 


6.85 


0,861 99 


0.855 96 


0,822 06 


0.901 44 


0,849 69 


0,873 70 


60 


16.85 


0,880 06 


0.874 84 


0.845 44 


0,915 01 


0.869 44 


0,890 52 


60 


36.85 


0,908 67 


0.904 46 


0.881 83 


0.936 74 


0,900 52 


0,917 23 


60 


56.85 


0,930 11 


0.926 96 


0,908 68 


0.953 1 8 


0.923 68 


0,937 30 


120 


-3.15 


0,721 33 


0.710 44 


0.641 45 


0.810 24 


0.695 40 


0,750 74 


120 


6.85 


0,760 25 


0,750 66 


0,689 71 


0.837 82 


0,737 80 


0,785 86 


120 


16.85 


0.793 17 


0,784 75 


0.731 23 


0.861 37 


0,773 69 


0,815 69 


120 


36,85 


0.845 15 


0,838 63 


0.796 98 


0.899 13 


0.830 22 


0,66311 


120 


56.85 


0.883 83 


0,878 70 


0.845 53 


0.927 66 


0,87211 


0.898 62 
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Annex D 

(normative) 

Pressure and temperature conversion factors 

If the input variables pressure and temperature are not in the necessary units of megapascals and kelvins, then 
conversions will have to be made in order to use the Fortran implementation. A selection of appropriate conversion 
factors is given below: 

Pressure 

p(MPa) = {p{bar)lxlO-i 

/> (MPa) = tp(atm)) X 0.101 325 

p (MPa) » tp(psia)yi 45,038 

p (MPa) s (A)(psig) + 14,695 9]/145,038 

Temperature 

r(K) = /CO) + 273,15 
T (K) = t'CP) - 32)/1 ;a + 273,15 
r(K) = lt("R)l/1.8 



16 



©ISO 



ISO 1221 3-2:1 997(E) 



Annex E 

(informative) 

Performance over wider ranges of application 

The AGA8-92DC equation has been comprehensively tested, over the temperature range 263 K to 338 K and at 
pressures up to 30 MPa, with the GERG data PI, (31 and the Gas Research institute data W for gases within the 
ranges of composition given for pipeline quality gases (see 4.4.1). Within these limits, the uncertainties are as given 
in 4.5. 

Rough estimates of the uncertainties involved in calculations of compression factors for wider ranges of application 
(with respect to composition) (see 4.4.2) are plotted in figures E.1 to E.4 as pressure-composition plots for nitrogen, 
carbon dioxide, ethane and propane, respectively. 

• In figures E.1 to E.4 the perfomiance of the AGA8-920C method is illustrated up to a maximum pressure of 
30 MPa. The uncertainty limits are dependent upon pressure, temperature and composition, and are also strongly 
affected by the proximity of the phase boundary. The estimated uncertainty limits presented below are based upon 
less comprehensive data, published as a supplement to the GERG databank (3], and upon the databanks in 
references [2] and [4], Reference [4] also provides data up to 70 MPa. The uncertainty limits given in figures E.1 to 
E,4 are always for the worst-case result, l.e. they are the least optimistic choice. 

Dashed lines are used to separate two regions of estimated uncertainty when the experimental evidence is not 
sufficient to determine the position of the boundary. The detailed composition of the gas will have a strong influence 
on the position of the phase boundary and the user should, therefore, make his own phase boundary calculation. 
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AGA8-0C92 equation 




r-2&3Kto338K 


B 






A 


A A/ a t0.1 % 




B «0.1 % to % 
1 


1 '. 1 \ 



0 0.1 0.2 0.3 O.t O.S 

Mole fraction of nitrogen, x^j 

Figure E.l — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high nitrogen content 
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Figure E.2 — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high carbon dioxide content 
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B £0.1 % to £0.2 % 
C £0.2 % to £0.5 % 



0.05 



0.1 



0.1S 0.2 

Mole fraction of ethane. Xc?^ 



Figure E.3 — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high ethane content 
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A6A&-0C92 equation 
r«263Kto 338K 



A A^atO,1% 
B £0.1%(ox0.2% 



O.OS 



0.1 

Mole fraction of propane. Xfiu 



Figure E.4 — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high propane content 

The overall results at pressures up to 10MPa and temperatures within the range 263 K to 338 K can be 
summarized as follows. Only gsises having mole fractions within the limits given below will have uncertainties within 
±0,1 %, ± 0,2 % and ± 0,5 %, respectively, within the given pressure and temperature domain. 



Component 


Mole fraction for an uncertainty within 




±0.1 % 


±0.2% 


±0.5% 


Nitrogen 


«0.50 






Carbon dioxide 


<0,23 


<s0.26 


<0.28 


Ethane 


<0,13 


<0.20 




Propane 


<0.06 


^0,10 
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Annex F 

(informative) 

Subroutines in Fortran for the AGA8-92DC method 

c 

C update: 17.05.94 E.W. Lemmon/S.W. Beyerlein/J.L. Savidge 

C update: 4.09.95 M. Jaeschke J. Sikora 

C AGA8-DC92 COMPRESSION FACTOR EQUATION - 

C SUBROUTINE DCAGA 

C • - 

C This progr€ua was written to accompany ISO 12213. 
C 

C "DCAGA* Calculates the compression factor of natural gases using 
C a detailed gas analysis. 

C 

C For information contact: DR. Jeffrey L. Savidge 
C gas research institute 

C 8600 W. Bryn Mawr Ave. 

C Chicago, IL 60631 

C (312) 399-8100, FAX (312) 399-8125 

C 

^ C This program calculates compression factors and molar densities for 

g, C natural gases from the input of gas conposition in accordance with the 

g C AGA8-DC92 compression eq[uation developed by the Gas 

S C Research Institute, Chicago, Illinois. (K.E. Starling and J.L. 

C Savidge, Compressibility Factors of Natural Gas and Other Related 

g C Hydrocarbon Gases, American Gas Association, AGA Transmission 

« C Measurement Committee Report No. 8, American Petroleum Institute MPMS 

2 C Chapter 14.2, Second Edition, November 1992, Catalog No. XQ9212) . 

t c 

S C The coefficients used in this program are the same as the values found 

« C in AGA Report No. 8, November 1992. 

§ C Values for the gas constant and molar masses conform with ISO 6976 

^ C (1995) and GPA 2172 (1988). 

2 C ■ 

1 ^ 

< C Ranges of application for compression factor calculation 
C with the AGA8 - DC92 equation: 

C 

C p-T-ranges 

C absolute pressure 0 to 65 MPa 

C tenperature 225 to 350 K 

C 

C Ranges for percentage molar composition: • 

C 

C A: pipeline quality natural gas 

C B: wider ranges of application 

C 



C 






A 








B 






c 


methane 


70 


to 


100 




50 


to 


100 




c 


nitrogen 


0 


to 


20 




0 


to 


50 




C 


carbon dioxide 


0 


to 


20 




0 


to 


30 




C 


ethane 


0 


to 


10 




0 


to 


20 




C 


propane 


0 


to 


3. 


5 


0 


to 


5 




c 


butanes 


0 


to 


1. 


5 


0 


to 


1, 


5 


c 


pentanes 


0 


to 


0. 


5 


0 


to 


0. 


5 


c 


hexanes 


0 


to 


0. 


1 


0 


to 


0. 


1 


c 


heptanes > 


0 


to 


0. 


05 


0 


to 


0. 


05 


c 


octanes plus 


0 


to 


0, 


05 


0 


to 


0. 


05 


c 


hydrogen 


0 


to 


10 




0 


CO 


10 




c 


carbon monoxide 


0 


to 


3 




0 


to 


3 




c 


helium 


0 


to 


0. 


5 


0 


to 


0. 


5 
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C water 0 to 0.015 0 to 0,015 

C 

C The expected uncertainty o£ the calculated results are for 

C pipeline quality natural gases: 

C +/- 0,1% within p-T range 0 to 12 MPa, 263 to 350 K 

C +/- 0.6% within p-T range 0 to 60 MPa, 225 to 350 K 

C The expected xincertairity in the wider remge of application 

C (composition) is often even for pressures up to 12 MPa larger. 

C For more details see ISO 12213 part 2, 

C 

SUBROUTINE DCAGA (XJ) 

INTEGER B(58) ,C(58) ,K(58} ,0(58) 
INTEGER Q(58) ,F{58) ,S{58) ,W(58) 
REAL*8 A(58) ,U(58) 

COMMON /CONSTANTS/ A, B , C, K, U,G,Q,-F, S, W 

REAL* 8 MW(21) ,EI(21) , KI (21) , GI (21) ,QI{21) ,FI(21) ,SI(21) ,WI(21) 

REAL*8 EIJ(21,21) .UIJ(21,21) ,KIJ (21 . 21) , GIJ(21 , 21) 

COMMON /PARAMETERS/ MW, EI, KI,GI,QI, FI, SI , WI.EU, UIJ, KIJ,GIJ 

REAL*8 Kl, CNS{58), 31(18) 
COMMON /COEF/ Kl , CNS. BI 

REAL*8 MWX, RGAS, TCM. DCM 
COMMON /MW/ MWX, RGAS, TCM; DCM 

INTEGER I. J, N 
REAL*8 SUM, XI (21), XJ(21) 
REAL*8 Ul, Gl, Ql, Fl, El 
REAL*8 XIJ, EIJO, GIJO , BN 

XI( 1) = XJ( 1) 

XI( 4) = XJ( 2) 

XI( 5) = XJ( 3) 

XI(ll) = XJ( 4) 

XI (12) = XJ( 5) 

XI (13) = XJ( 6) 

XI{14) « XJ( 7) 

XI (15) = XJ( 8) 

XI(16) = XJ( 9) 

XI(17) = XJ(IO) 

XI(18) = XJ(ll) 

XI(19) = XJ(12) 

XI( 3) = XJ(13) 

XI( 2) = XJ(14) 

XI( 7) = XJ(15) 

XI(20) = XJ(16) 

XI( 6) = XJ(17) 

XI(IO) = XJ(18) 

XI(21) = XJ(19) 

XI( 8) = XJ(20) 

XI( 9) = XJ(21) 

C Normalize mole fractions 

SUM = 0 

MWX = 0 

DO 10 1=1, 21 

10 SUM = SUM + XI (I) 

DO 20 1=1, 21 

20 XI(I) = XI (I) /SUM 

C Calculate molecular weight 

RGAS = 8.31451D-3 
MWX = 0 
DO 30 1=1, 21 
30 MWX = MWX + XI(I)*MW(I) 

DO 40 N=l, 18 
40 BI(N) = 0 
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Kl = 0 

Ul = 0 

Gl = 0 

Ql = 0 

Fl := 0 

El = 0 



DO 50 1=1, 21 

Kl = Kl + XI(I)*KI{I)**2.5D0 
Ul = Ul XI(I)*EI{I)**2.5D0 
Gl = Gl 4- X1(X)*GICI) 
Ql = Ql + XI(I)*QI(I) 
Fl = Fl + XI(I)*XI(I}*FI(I) 
El = El XI(I)*EI(I) 
50 CONTINUE 

TCM - 1.261*E1 

DCM = K1**{-1.2D0) 

Kl = K1*K1 

Ul = U1*U1 



• DO 60 1=1, 8 

DO 60 J=^I-fl, 19 

' XIJ = Xi{I)*XI(J) 

IF (XIJ.NE-0) THEN 

Kl = K1+2.D0*XIJ*{KIJ{I,J)**5.D0-1,D0)*(KI(I)*KI(J))**2.5D0 

Ul = U1+2.D0*XIJ*{UIJ(I,J)**5.D0-1.D0)*{EI(I)*EI(J))**2.5D0 

3 Gl = G1+XIJ*(GIJ(I,J) - l.DO)*(GI(I) + GI(J)> 

I ENDIF 

a 60 CONTINUE 
c 

I DO 80 1=1, 21 

3 DO 80 J=I, 21 

e XIJ = XI(I)*XI(J) 

IF (XIJ.NE.O) THEN 

IF (I.NE.J) XIJ = 2.D0*XIJ 
H EIJO = EIJ(I,J)*DSQRT(EI{I)*EI(J)) 

^ GIJO = GlJd, J) * (GI(I) + GI(J))/2.D0 

^ DO 70 N=l, 18 

S BN = (GIJO 4. l.DO - G{N))**G{N) 

§ & * (QI(I)*QI(J) + l.DO " Q(N))**Q(N) 

^ & * (DSQRT(FI (I) *FI(J)) f l.DO - F(N))**F(N) 

H & * (SI(I)*SI(J) + l.DO - S(N))**S(N) 

« & * (WI{I)*WI{J) + l.DO - W{N))**W(N) 

< BI(N) = BI(N)+A(N)*XIJ*EIJO**U{N) *(KI{I)*KI{J) )**1.5D0*BN 

70 CONTINUE 

ENDIF 
80 CONTINUE 



•9 



Kl = K1**0.2DO 

Ul = U1**0.2DO 

DO 90 N=13, 58 
90 CNS(N) = (Gl + l.DO - G{N))**G(N) 

& * (Ql**2 + l.DO - Q(N))**Q(N) 

& * (Fl l.DO - F{N))**F(N) 

fie * A(N) *U1**U(N) 

END 

SUBROUTINE P20FDT(D, T, P, Z, BMIX) 

INTEGER B(58) ,C(58) ,K{58) ,G{58) 
INTEGER Q(58) ,F(58) ,S{58) ,W(58) 
REAL*8 A{58) ,U(58) 

COMMON /CONSTANTS/ A, B. C, K, U, G,Q, F. S, W 

REAL*8 Kl, CNS(58), BI(18) 
COMMON /COEF/ Kl, CNS , BI 

REAL*8 MWX, RGAS, TCM, DCM 
COMMON /MW/ MWX, RGAS, TCM, DCM 
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" INTEGER N 
REAL*8 T. P, Z, BMIX. DR 

DR = D*K1**3 
BMIX =0 
DO 10 N=l, 18 
10 BMIX- « BMIX BI(N)/T**U(N) 

2 = l.DO + BMIX*D 
DO 20 N=13, 18 
20 Z = Z - DR*CNS(N)/T**U(N) 



DO 30 N=13, 58 

30 Z = 2 + CNS(N) /T**U(N) * (B{N) - C (N) *K (N) »DR**K(N) ) ♦DR**B (N) 
& •DEXP(-C{N)*DR**K(N) ) 

P = D*RGAS*T*2 
END. 

SUBROUTINE D20FPT(P, T, D, Z, BMIX) 

REAL*a P, T, D, 2, BMIX 

REAL*8 XI, X2. X3 , F, Fl. F2 , F3, TOL 

INTEGER I 

TOL = 0.5D-9 
XI = O.OOOOOIDO 
X2 = 40. DO 
D= 0 

CALL P20FDT(X1. T, Fl, 2, BMIX) 
CALL P20FDT{X2, T, F2 , 2. BMIX) 
Fl = Fl - P 
F2 = F2 - P 

IF (Fl*F2.GE.O) RETURN 



C 
C- 



BEGIN ITERATING 



DO 60 I = 1, 50 

...Use False Position to get point 3. 
X3 = XI - F1*(X2 - XI) /(F2 - Fl) 
CALL P20FDT(X3, T, F3 . 2, BMIX) 
F3 = F3 - P 



C ...Use points 1, 2, and 3 to estimate 
C . . .method (quadratic solution) . 



the root using Chamber's 



Xl*F2*F3/( (Fl 
& + X2*Fl*F3/( (F2 
& + X3*Fl*F2/( (F3 
IF ( (D - XI) * (D - 
CALL P20FDT(D, 
F = F - P 

IF (DABS (F) .LE. TOL) RETURN 



F2) *(F1 

- F1)*(F2 

- F1)*(F3 
X2) .GE.O) D 
F, Z, BMIX) 



F3)) 
F3)) 
F2)) 
= (XI 



X2) /2.D0 



C ...Discard quadratic solution if false position root is closer. 
IF (DABS (F3 ) .LT. DABS (F) .AND. F*F3.GT.O) THEN 
IF (F3*F1.GT.0) THEN 
XI = X3 
Fl = F3 
ELSE 

X2 = X3 
F2 = F3 
ENDIF 
ELSE 



..Swap in new value from quadratic solution 
IF (F*F3.LT.O) THEN 
XI = D 
Fl = F 
X2 = X3 
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F2 = F3 
ELSEIF (F3*F1.GT.0) THEN 
XI = D- 
Fl •= F 
ELSE 
X2 = D 
F2 = F 

ENDIF . • 

■ ENDIF 
60 CONTINUE 
D = 0 
END 

BLOCK DATA 

INTEGER B(58),C(58),K(58),G(58) 
INTEGER Q(58) ,F(58y,S{58) ,W(58) 
REAL*8 A(58) ,U{58) 

COMMON /CONSTANTS/ A, B, C, K,U,G. Q, F,S,W 

REAL*8 MW(21) ,EI(21) .KI(21) .GI(21) .QI{21) ,FI(21) ,SI(21) ,WX(21) 

REAL*8.EIJ(21,21) ,UIJ(21,21) ,KIJ(21,21) ,GIJ{21,21) 

COMMON /PARAMETERS/ MW, EI , KI, GX,QI , FI , SI , WI , EIJ.UIJ, KIJ ,GIJ 



REAL*8 XN{21), XH(21) 
COMMON /GRENZDATA/ XN . XH 



- C Equation of state parameters 

g DATA A/ 

& 0,153832600D0, 
& 0.375796500DO, 
& -0,71023704000, 
& 0.2291290D-08, 
& -0.003433888DO, 
fic -0.001600573D0, 
& 0.279179500D0, 
& 3 .150547000DO, 
& -0.079314910DO, 
& 0.034317290D0, 
& 0.746545300DO, 
& -0.005399808DO, 
& X.874951000DO, 
& 0.009776195D0, 
& -0.001226752D0, 



tit 
c 

3 
I 



1.341953000D0. -2.99858300000, 
-1.589575000DO, -0 . 053588470DO , 
-1.471722000DO. 1 . 321850350DO , 
0.157672400DO, -0 . 4363 86400D0, 
0.032059050DO, 0 . 024873550DO , 
0.642470600D0, -0 . 416260100DO , 
-0.696605100D0, -0 . 002860589D0 , 
0.007224479DO, -0 . 70575290OD0, 
-1.418465000DO. -0 . 5999050D-16, 



-0.O48312280D0 / 
0.886594630DO, 
-0,78665925000, 
-0.044081590DO, 
0.073322790D0, 
-0.066899570DO, 
-0.OO8O98836D0, 
0.534979200D0, 
0.105840200D0, 
0.042968180DO, 
-9.936757000D0, 
0.003733797D0, 
0,000205518DO, 
0,00686241500, 



-0. 007022 847D0, 0 . 024955870DO , 
-0.291961300DO, 7 . 294616000D0 , 
-0,24325670000. 0 . 049870160D0 , 
0.002168144D0, -0 . 658716400DO , 
-0.020487080D0, 0 . 015573220DO, 
0,00285090800/ 

DATA B/1, 1.1, 1,1, 1,1. 1,1, 1,1, 1,1. 1,1, 1,1, 1,2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 
& 3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,5,5,5,5,5,6,6,7,7,8,8,8,9,9/ 

DATA C/0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 1,1, 1,1, 1,0, 0,1, 1,1, 1,1, 1,1,0,1, 
& 1,1,1,1,1,1,1,1,0,0,1,1,1,1,1,0,1,1,1,1,0,1,0,1,1,1,1,1,1/ 

DATA K/0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,3, 2, 2, 2. 4, 4, 0,0, 2, 2, 2, 4,4,4, 4,0,1, 
& 1,2,2,3,3,4,4,4,0,0,2,2,2,4,4,0,2,2,4,4,0,2,0,2,1,2,2,2,2/ 

DATA G/ 0,0, 0,0, 1,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,1, 0,0, 0,1, 
& 0,0, 1,1, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0. 0,0, 0,1, 0,0, 1,0, 1,0,0/ 

DATA Q/0,0, 0,0,0,0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,1,0,1,0, 
& 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 1,0. 0,0, 0,1, 0,1, 0,0, 1,0, 0,0, 0,0,1/ 

DATA F/0. 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 
& 1,0, 0.0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0.0, 0,0, 0,0, 0,0, 0,0,0/ 

DATA S/0, 0,0, 0,0, 0,0, 1,1. 0,0. 0.0, 0,0. 0.0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
& 0,0, 0,0, 0.0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0.0. 0,0, 0,0, 0.0, 0,0,0/ 

DATA W/0, 0,0, 0.0, 0,0, 0,0, 1,1, 1,0, 0,0, 0.0. 0.0, 0,0, 0,0. 0,0, 0,0, 0,0, 
& 0,0, 0,0,0,0,0, 0,0, 0,0, 0,0. 0.0. 0.0. 0,0. 0.0. 0.0. 0,0, 0.0. 0,0/ 

DATA U/0 , 0 . 5D0 , 1 , 3 . 500 , - 0 . 500 , 4 . 5D0 , 0 , 5D0 , 7 . 5D0 , 9 . 500 , 6 , 12 . 12 . 5D0 , 
& -6.2.3,2,2,11.-0.500,0.500,0,4,6,21,23,22,-1,-0.500,7,-1,6, 
& 4, 1,9, -13,21,8,-0.500,0.2.7,9,22,23. 1.9.3,8,23,1.500,5, 

& -0.500,4,7.3,0,1,0/ 



, .Characterization Parcixneters 

DATA MW/16. 043000, 28.0135D0, 44.010000, 30.070000, 44.097000, 

& 18.015300, 34.082000, 2. 015900, 28.0100DO, 31.998800, 

t 58.1230D0, 58,123000, 72.1500DO, 72.150000, 86.177000, 

& 100.204000,114,231000,128.258000.142.285000, 4.0026DO. 
& 39.948000/ 
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DATA EI/151. 318300D0. 99.737780D0, 241 . 960600DO , 244 . 166700D0 . 

& 298.118300D0, 5 14 . 015600D0 . 296 .355000D0, 26.957940D0, 

& • 105.534800D0, 122 . 766700D0, 324 . 068900DO , 337 . 638900D0 , 

&. 365.599900D0. 370.682300DO, 402 . 636293D0 , 427 .722630D0 , 

& 450,32502200, 470 . 840891D0, 489 . 558373D0 , 2.610111D0, 

& 119.629900D0/ 

DATA KI/0.4619255D0, 0.4479153D0, 0.4557489D0, 0.5279209D0. 

& 0.5837490D0, 0.3825868D0, 0.4618263D0, 0.3514916D0, 

& 0.4533894D0, O.4186954D0, 0.6406937D0, 0.6341423D0. 

& 0.673857700/ 0.6798307D0, 0. 717511800, 0. 752518900, 

& 0.7849550D0, 0. 815273100, 0. 843782600. 0. 358988800, 

& 0.421655100/ 

DATA GI/0, 0.027815D0, 0.189O65D0, 0.079300D0, 0.141239D0, 

& 0. 33250000, 0. 08850000, 0. 03436900, 0,03895300, 0.021000DO, 

& 0. 25669200, 0. 28183500, 0.332267O0, 0. 36691100, 0.289731D0, 

& 0. 33754200, 0. 38338100, 0. 42735400, 0,46965900, 0, 0/ 

DATA QI/2*0, 0.69O0, 2*0, 1.06775D0, O.633276D0, 14*0/ 



DATA FI/7*0, 1, 13*0/ 

DATA SI/5*0, 1,582200, 0.390O0, 14*0/ 
DATA WI/5*0, 1, 15*0/ 

C Binary interaction parameters 

DATA EIJ/441*1/ 
DATA UIJ/441*1/ 
DATA KIJ/441*!/ 
DATA GIJ/441*1/ 



DATA (EIJ(1, J) , J=2,19) / 
& 0. 97164000, 0,96064400, 1, 0.994635D0, 0.708218DO, 

& 0.931484DO, 1.170520D0, 0.990126D0, 1, 1,01953000, 

& 0.98984400, 1. 00235000, 0,99926800, 1. 10727400, 0. 88088000, 

& 0. 88097300, 0,88106700, 0.881161D0/ 

DATA (EIJ(2, J) , J=3,141 / 
& 1.022740D0, 0,97012000, 0. 94593900, 0.74695400, 0. 90227100, 

& 1.086320DO, 1,00571000, 1. 02100000, 0. 94691400, 0. 97338400, 

& 0. 95934000, 0.94552000/ 

DATA (EIJ(3, J) , J=4,19) / 
& 0.925O53D0, 0. 96023700, 0. 84940800, 0. 95505200, 1.281790D0, 

& 1.5O0, 1, 0,90684900, 0.897362O0, 0.72625SO0, 

& 0. 85976400, 0.855134D0, 0.83122900^ 0. 80831000, 0.786323D0, 

& 0,76517100/ 

DATA (EIJ(4, J) ,J=5,14) /1.022560O0, 0. 69316800, 0. 94687100, 
& 1.164460D0, 3*1, 1.O13O60D0, 1 1 , 0053200/ 

DATA (EIJ(5, J) ,J=8,12)/i. 03478700, 3*1, 1.004900/ 

DATA (EIJ(7,J) ,J=15,19)/1.008692DO, 1,01012600, I.OIISOIDO, 
& 1. 01282100, 1.014089D0/ 

DATA {EIJ(8, J) , J=9,12)/1.1D0, 1, 1.3O0, 1.300/ 

DATA (UIJ(1, J) , J=2,19) / 

& 0.886106D0, 0. 96382700, 1, 0.990877DO, 1, 0. 73683300, 

& 1. 15639000, 3*1, 0. 99229100, 1, 1. 00367000, 1.30257600, 

& 1.191904D0, 1.205769D0, 1.21963400, 1.23349800/ 

DATA (UIJ(2,J) ,J=3,12)/O-835O58D0, 0. 81643100, 0.915502D0, 1, 

& 0. 99347600, 0. 40883800, 3*1, 0.99355600/ 

DATA (UIJ(3, J) ,J=4,19)/0. 96987000, 2*1, 1.045290D0, 1. 0 . 9O0 , 

& 5*1, 1-066638DO, 1.077634D0, 1.088178D0, 1.098291D0, 

& 1.108021D0/ 

DATA (UIJ{4, J) , J=5,14) / 

& 1. 06517300, 1, 0.971926D0, 1.61666000, 2*1, 4*1,2500/ 

DATA {UIJ(7, J) , J=7,19)/8*l, 1. 02897300, 1.033754DO, 

& 1. 03833800, 1. 04273500, 1.04696600/ 

DATA (KlJd, J) ,J=2,19)/ 
& 1.003630D0, 0. 99593300, 1, 1.007619DO, 1, 1. 00008000, 

& 1,02326000, 3*1, 0,99759600, 1, 1.00252900, 0.98296200, 

& 0.983565D0, 0,98270700, 0. 98184900, 0.98099100/ 

DATA (KIJ(2,J) ,J=3,8)/ 
& 0.982361D0, 1.007960DO, 1, 1, 0.942596000,1.03227000/ 

DATA (KIJ{3,J) .J=4,19)/1,003510D0, 2*1, 1.0077900,7*1.000, 
& 0.910183DO, 0. 89536200, 0.881152D0. 0.867520DO, 0.85440600/ 
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T3 
C 

E 
tJ 

a 
c 



c 

3 

U 
3 

I! 

3 
< 



DATA (KIJ(4,J),J=>5,8)/0.986893D0, 1, 0.999969D0, 1.020340D0/ 
DATA (KIJ(7,J) ,J=7,21)/8*1, 0.968130D0, 0.962870D0, 
& 0.957828D0, 0.952441D0, 0.948338D0, 2*1/ 

DATA GIJ(1,3) /0.807653D0/ 
DATA GIJ(1,8). /1.957310D0/ 
DATA GIJ{2,3) /0.982746DO/ 
DATA GIJ{3,4) /0.370296D0/ 
DATAGIJ(3,6) /1.673090D0/ 



DATA XN/ 
& 50.0 . 
& 0.0 . 
& 0.0 , 



0.0 
0-0 
0.0 



0.0 
0.0 
0.0 



0-0 
0,0 
0.0 



0.0 
0.0 
0.0/ 



0.0 
0.0 



0-0 
0.0 



O.Q 
0.0. 



DATA XH/ 
& 100.0 , 
& 0.05 , 
& 0.015, 

END 



20.0 
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0.02 



5.0 

0.05 

0,02 



1.5 
0.05 
10.0 



. 1.5 
30.0 
3.0/ 



0.5 
50.0 



0.5 
0.02 



0.1 
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Foreword 

ISO (the International Organization for Standardization) is a worldwide 
federation of national -standards bodies (ISO member bodies). The work of 
prepanng Internationa! Standards is normally carried out through ISO 
technical committees. Each member body interested in a subject for which 
a technical committee has been established has the right to be represented 
on that committee. International organizations, govemmental and non- 
governmental, in liaison with ISO, also take part In the work. ISO 
collaborates closely with the International Electrotechnical Commission 
(lEC) on all matters of electrotechnical standardization. 

Draft Intemational Standards adopted by the technical committees are 
circulated to the member bodies for voting. Publication as an Intemational 
Standard requires approval by at least 75 % of the member bodies casting 
a vote. 

International Standard ISO 12213-3 was prepared by Technical Committee 
ISO/TC 193, Natural gas. Subcommittee SC 1, Analysis of natural gas, 

ISO 12213 consists of the following parts, under the general title Natural 
gas — Calculation of compression factor. 

— Part 1: Introduction and guidelines 

— Part 2: Calculation using molar-composition analysis 

— Part 3: Calculation using physical properties 

Annexes A to D form an integral part of this part of ISO 12213. Annexes E 
to G are for information only. 
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Natural gas — Calculation of compression factor — 
Part 3: 

Calculation using physical properties 

1 Scope 

This International Standard specifies methods for the calculation of compression factors of natural gases, natural 
gases containing a synthetic admixture and similar mixtures at conditions under which the mixture can exist only as 
a gas. 

This part of ISO 12213 specifies a method for the calculation of compression factors when the superior calorific 
value, relative density and carbon dioxide content are known, together with the relevant pressures and 
temperatures. If hydrogen is present, as is often the case for gases with a synthetic admixture, the hydrogen 
content also needs to be known. 

NOTE — In principle, it is possible to calculate the compression factor when any three of the parameters superior calorific 
value, relative density, carbon dioxide content (the usual three) and nitrogen content are known, but subsets including nitrogen 
content are not recommended. 

The method is primarily applicable to pipeline quality gases within the ranges of pressure p and temperature T at 
which transmission and distribution operations normally take place, with an uncertainty of about ± 0,1 %. For wider- 
ranging applications the uncertainty of the results increases (see annex E). 

More detail concerning the scope and field of application of the method is given in part 1 of this International 
Standard. 



2 Normative references 

The following standards contain provisions which, through reference in this text, constitute provisions of this part of 
ISO 12213. At the time of publication, the editions indicated were valid. All standards are subject to revision, and 
parties to agreements based on this part of ISO 12213 are encouraged to investigate the possibility of applying the 
most recent editions of the standards indicated below. Members of lEC and ISO maintain registers of currently valid 
International Standards. 

ISO 31-3:1 992, Quantities and units — Part 3: Mectmnics. 
ISO 31 -4:1 992. Quantities and units — Pari 4: Heat, 

ISO 6976:1995, Natural gas — Calculation of calorific values, density, relative density and Wobbe index from 
composition. 

ISO 1 221 3-1 :1 997, Natural gas — Calculation of compression factor — Part 1: Introduction and guidelines 
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3 Definitions 

All definitions relevant to the use of this part of ISO 12213 are given In part 1. 



4 Method of calculation 
4.1 Principle 

The method recommended uses equations which are based on the concept that pipeline quality natural gas may be 
-uniquely characterized for calculation of its volumetric properties by an appropriate and distinctive set of measurable 
physical properties. These characteristics, together with the pressure and temperature, are used as input data for 
the method. 

The method uses the following physical properties: superior calorific value, relative density and carbon dioxide 
content. The method is particulariy useful in the common situation where a complete molar composition is not 
available, but may also be preferred for its relative simplicity. For gases with a synthetic admixture, the hydrogen 
content needs to be known. 



4.2 The SGERG-88 equation 

The calculation method using physical properties is based on the standard GERG 88 (SGERG-88) virial equation 
for natural gases ["*!. [2]. {31. The standard GERG 88 virial equation is derived from the master GERG 88 
(MGERG-88) virial equation, which is a method of calculation based on a molar-composition analysts W, 

The SGERG-88 virial equation from which the compression factor Z is calculated may be written as 

Z = 1 + 5A„ + Cp^ ...(1) 
where 

B and C are functions of the input data comprising the superior calorific value Hq, the relative density d, the 
contents of tx)th inert and combustible non-hydrocarbon components of the gas mbcture (CO2 and H2) 
and the temperature T, 

Pm is the molar density given by 

p^^pjiZFT) ...(2) 

where 

^ =/i (p. 7; //s, d, XCO2. • • • (3) 

However, the SGERG-88 method treats the natural-gas mixture internally as a five-component mixture consisting of 
an equivalent hydrocarbon gas (with the same thermodynamic properties as the sum of the hydrocarbons present), 
nitrogen, carbon dioxide, hydrogen and carbon monoxide. To characterize the thermodynamic properties of the 
hydrocarbon gas adequately, the hydrocarbon heating value H^h i^ ^^^^ needed. Therefore, the calculation of Z 
uses 

In order to be able to model coke oven gas mixtures, the mole fraction of carbon monoxide is taken to have a fixed 
relation to the hydrogen content. If hydrogen is not present (jch2 < 0,001), then set = 0. The natural-gas mixture 
is then treated in the calculation method as a three-component mixture (see annex B). 

The calculation is performed in three steps. 

First, the five-component composition from which both the known superior calorific value and the known relative 
density can be calculated satisfactorily may be found from the input data by an iterative procedure described in 
detail in annex B. 
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Secondly, once this composition is known, B and C may be found using relationships also given in annex 6. 
In the third step, equations (1) and (2) are solved simultaneously for and Z by a suitable numerical method. 
A flow diagram of the procedure for calculating Z from the input data is shown in figure B.1 . 

4^ Input variables 

4.3.1 Preferred input data set 

The input variables required for use with the SGERG-88 equation are the absolute pressure, temperature and 
superior calorific value (volumetric basis), the relative density, the carbon dioxide content and the hydrogen content. 
Thus the physical properties used the input data set (set A) are 

^s. ^. Jfcoa and JCH2 

Relative density is referred to normal conditions (101,325 kPa and 0 ""C) and superior caloriflc value is referred to 
normal conditions (1 01 ,325 kPa and 0 "C) and a combustion temperature of 25 *C. 

4.3.2 Altemative Input data sets 

Three alternatives to the preferred input data set (see 4.3.1) may be used with the standard GERG virial equation: 
JCN2. Ws.*' and JCH2 (sets) 
- JTNa. Jfcoa- ^ and JCH2 (set C) 
XNzi J^coz. and xy^^ (set D) 

The altemative input data sets are considered fully in GERG Technical Monograph TM5 PI. Use of the alternative 
input data sets gives results which may differ at the fourth decimal place. This part of ISO 12213 recommends the 
use of input data set A. 

4.4 Ranges of application 
4.4.1 Pipeline quality gas 

The ranges of application for pipeline quality gas are as defined below: 

absolute pressure 0 MPa is 12 MPa 

temperature 263 K ^338K 

mole fraction of carbon dioxide 0 ^ jcco2 ^-^^ 

mole fraction of hydrogen 0 < x^j ^ 0,10 

superior calorific value 30 MJ nrrS =^45 MJ m-3 

relative density 0,55 < 0,60 

The nnole fractions of other natural-gas components are not required as input These mole fractions shall, however, 
fie within the following ranges: 



methane 


0,7 




1,0 


nitrogen 


0 




:SO,20 


ethane 


0 






propane 


0 




^ 0,035 


butanes 


0 




^ 0.015 
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pentanes 

hexanes 

heptanes 



octanes plus higher hydrocarbons 
carbon monoxide ' 



helium 



0 



0 



0 



0 



0 



0 



< xce ^ 0,001 

*s XC7 ^ 0,000 5 

"Sxce^ ^0,000 5 
^ xco 0,03 

«JCHe ^0,005 



water 



0 



^JrH20 0,000 15 



The method applies only to mixtures in the single-phase gaseous state (above the dew point) at the conditions of 
temperature and pressure of interest For pipeline quality, the method is applicable over wider ranges of 
temperature and pressure but with increased uncertainty (see figure 1). In the computer implementation, the lower 
temperature limit is set at 250 K. 

4.4.2 Wider ranges of application 

The ranges of application tested beyond the limits given in 4.4.1 are: 

absolute pressure 0 MPa ^12 MPa 

temperature 263 K =s338K 

mole fraction of carbon dioxide 0 ^ xcoz ^ 0,30 

mole fraction of hydrogen 0 < jch2 < 0,10 

superior calorific value 20 M J-nrr3 ^Hs ^48 MJ nrr^ 

relative density 0,55 ^ 0,90 

The allowable mole fractions of other major natural-gas components are extended to: 

methane 0,5 ^arou ^1,0 

nitrogen 0 ^ 0,50 

ethane 0 ^ XC2H5 0,20 

propane 0 ^ xcaHa 0,05 

The limits for other minor natural*gas components remain as given in 4.4.1 for pipeline quality gas. 

The method is not applicable outside these ranges; the computer implementation described in annex 6 will not allow 
violation of the limits of composition quoted here. 



4^ Uncertainty 

4,5.1 Uncertainty for pipeline quality gas 

The uncertainty in the prediction of the compression factor AZ (for the temperature range 263 K to 338 K) is ± 0.1 % 
at pressures up to 10 MPa and ±0,2% l>etween 10 MPa and 12 MPa for natural gases with xn2^0,20. 
xco2^0,09, xc2Hb^0;10 and xh2^0,10, and for 30 MJ-nr^ ^ Hs 45 MJ-m-3 and 0,55 0,80. (see 
figure 1). 

For gases with a CO2 content exceeding 0,09, the uncertainty of ± 0,1 % is maintained for pressures up to 6 MPa 
and for temperatures between 263 K and 338 K. This uncertainty level is determined by comparison with the GERG 
databank on measurements of the compression factor for natural gases l^l. 161 and with the Gas Research Institute 
data [9). 
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350 

Temperature, r (K) 



Figure 1 — Uncertainty limits for the caiculation of compression factors (The uncertainty limits given are 
expected to be valid for natural gases and similar gases with < 0,20; jccq2 ^ 0.09; jrc2H6 ^ 0,10 and jch2 < 0,10, 

and for 30 MJrn-3 //s 45 MJnr3 and 0,55 0,80) 

4.5.2 Uncertainty for wider ranges of application 

The estimated uncertainties involved in calculations of compression factors beyond the limits of quality given in 
4.5.1 are discussed in annex E. 

4.5.3 Impact of uncertainties of input variables 

Listed in table 1 are typical values for the uncertainties of the relevant input variables. These values may be 
achieved under optimum operating conditions. 

As a general guideline only, an error propagation analysis using the above uncertainties in the input variables 
produces an additional uncertainty of about ± 0,1 % in the result at 6 MPa and within the temperature range 263 K 
to 336 K. Above 6 MPa, the additional uncertainties are greater and increase roughly in direct proportion to the 
pressure (see reference [3])^ 



4.5.4 Reporting of results 

Results for the compression factor shall be reported to four places of decimals, together with the pressure and 
temperature values and the calculation method used (ISO 12213-3, SGERG 88 equation). For verification of 
calculation procedures, it is useful to carry extra digits. 
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Table 1 — Uncertainties of input variables 



Input variable 


Absolute uncertainty 


Absolute pressure 
Temperature 

Mole fraction of cartx>n dioxide 
Mole fraction of hydrogen 
Relative density 
Superior calorific value 


± 0.02 MPa 
±0,15 K 
±0,002 
±0,005 
±0,001 3 
±0.06 MJm-3 



5 Suppliers of computer programmes 

it is planned to make software available whfcii implements this International Standard. Users are invited to contact 
their ISO member body or ISO Central Secretariat to enquire about the availability of such software. 
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Annex A 

(normative) 

Symbols and units 

The symbols specified in this annex are those which are used in both the main text and in annex B. The units 
specified here are those which give consistency with the values of the coefficients given in annex B. 



Symbol 

ftwo(O) 

bij(0) 

bij(2) 
B 

<^m 

CH2 

c«o(0) 

cm>(2) 
Cffi(0) 

c//i(2) 

C«2(1) 
c//2(2) 

Cy*(0) 
Cijk(2) 

C 



Meaning 

Zero-order (constant) term In the molar heating value (Wch) expansion of Bn 
[equation (B^O)] 

First-order (linear) term in the molar heating value (Hch) expansion of 
[equation (B^O)] 

Second-order (quadratic) term In the molar heating value (Hch) expansion 
of ^11 [equation (B.20)] 

Terms in the temperature expansion of bno [equation (B.21)] 
Terms in the temperature expansion of tf/i [equation (B.21)J 
Tenns in the temperature expansion of b^z [equation (B.21)l 

Terms in the temperature expansion of [equation (B.22)] 
Second virtal coefficient [equation (1)] 

Second virial coefficient for binary interaction between component i and 
component / [equation (B.22)] 

Zero-order (constant) term in the molar heating value (//ch) expansion of Cm 
[equation (B.29)] 

Hrst-order (linear) term in the molar heating value {Hch) expansion of C-m 
[equation (B.29)] 

Second-order (quadratic) term in the moiar heating value {Hch) expansion 
of Cm [equation (B.29)] 

► Terms in the temperature expansion of [equation (B.30)] 
I Terms in the temperature expansion of c//i [equation (B-30)] 
I Terms in the temperature expansion of [equation (8.30)] 

I Terms in the temperature expansion of Cy [equation (B.31)] 
Third virial coefficient [equation (1 )] 



Units 

m^kmohi 

m3.MJ-i 
m3.kmolMJ-2 

m3.kmol-i 

m3.kmoKi.K-i 

m3.kmol-iK-2 

m3.MJ-i 

m3.MJ-i-K-i 

m3-MJ-l.K-2 

m3.kmol-MJ-2 

m3.kmolMJ-2.K-i 

m3.kmol-MJ-2.K-2 

m3.kmohi 

m3.kmoH.K-i 

m3.kmoH.K-2 

m3.kmoH 
m3.kmoH 

m6.kmol-2 

mCkmoHMJ-t 

m6.MJ-2 

m6-kmoh2 

m6-kmol-2-K-i 

m6.kmol-2.K-2 

m6.kmoH-MJ-i 

m6kmoH.MJ-iK-i 

m6.kmoH.MJ-iK-2 

m6.MJ-2 

m6.MJ-2.K-i 

m6.MJ-2.K-2 

m6.kmol-2 

m6.kmoh2.K-i 

m6.kmoh2-K-2 

m5.kmol-2 
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Symbol . Meaning Units 

Cyjk Third virial coefficient tor ternary interaction between components /,y and Jt m^ kmoh^ 

[equation (B.31)l 

d Relative density [£/(air) = 1; equation (B.I )1 — 

DHcH Change in the molar heating value Hch during iteration ' MJ-kmoI-i 
[equations (B.10) and (B.11)] 

Hs Superior calorific value (gas at normal conditions (0 *'C» 1,013 25 bar), MJ nr3 

combustion temperature 25 '^C] 

H Molar heating value (combustion temperature 25 **C) . MJ kmol-i 

M Molar mass [equations (B.5) and {B.8)J kg kmohi 

p Absolute pressure bar 

R (Universal) gas constant rnS-bar-kmoH-K-i 

T Absolute temperature K 

/ Celsius temperature [=7-273.15; equation (B27)] **C 

Vjn Molar volume (= 1/pn^) rT>* kmol-i 

X Mole fraction of a component — 

y Combination rule parameters for the binary unlike-interaction virial coefficients — 

Biz ^nd Bi3 (table B.1) and the temary unlike-interaction virial coefficient C^^ 

[equation (B.32)] 

2 Compression factor — 

p Mass density [equations (B.8) and (B.42)l kg m-3 

Molar density (= Vin"^) kmol-mr^ 

Additional subscripts 

n Value at nomnal conditions (7n = 273,1 5 K, pn = 1 ,01 3 25 bar) 

CH For the equivalent hydrocarbon 

CO For carbon monoxide 

CO2 For cartDon dioxide 

H2 For hydrogen 

N2 For nitrogen 

Additional qualifiers 

(air) For dry air of standard composition [equation (B.I)] 

(D) For special value of p used in equation (B.11) 

1 For the equivalent hydrocarbon [equations (B.1 2) and (B.I 5)] 

2 For nitrogen [equations (B.I 2) and (B.I 6)] 

3 For carbon dioxide [equations (B.12) and (B.17)} 

4 For hydrogen [equations (B.I 2) and (B.I 8)] 

5 For carbon monoxide [equations (B.12) and (B,19)J 
(id) Ideal gas state 

(u) Iteration counter (B.2.1) 

(v) Iteration counter (B.2.2) 

(w) Iteration counter (B.4) 
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Annex B 

(normative) 

Description of the SGERG-88 method 

This annex gives the equations lor and numerical values of coefficients which together specify completely the 
SGERG method for calculation of compression factors. 

It also describes iteration procedures adopted by GERG (^l for implementing the method in the verified Fortran 77 
subroutine SGERG.FOR. This subroutine provides the correct solution; other computational procedures are 
acceptable provided that they can be demonstrated to yield identical numerical results. The calculated results shall 
agree to at least the fourth place of decimals with the examples given in annex C. 

Other implementations which are known to produce identical results are as follows: 

a) A Basic version, described in GERG TM5 PI, which may be used with a variety of metric reference conditions. 
This programme was designed mainly for PC applications. 

b) A version in C. described in German DVGW Directives, sheet G486 

c) A version in Turbo Pascal. 

All these programmes have been verified to give the same results to within 10-^. The availability of the programmes 
and the conditions which apply to their use are discussed in part 1 of this International Standard. 



B.I Basic structure of the calculation method 

As described in 4.2, the calculation proceeds in three steps, which are shown schematically in figure B.1. 



Input data: 










Calculate intermediate data: r„ • 273.1S K 
Pft. ^co. xoiiul XhjIoI. Hoiiu) 

(see clause 8.2) 










Calculate viriai coefficients: r « (/ ♦ 273,15) K 

(see clause B.3) 










Solve viriai equation: 

(see clause 8.4) 







Figure B.I — Flow diagram for standard GERG-88 calculation method 

(jt, = mole fraction of component i) 
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Input datat 




1 


Calculatei 


Pa. 


'CO 

Equations (B.D, (B.2) 






(/ * 0. 


V » 0 


1 




• Web » 0) 


-1000 . 
«-0.06S 






Calculate: 


Equations (8.3). (B.4) 












u = 




I 


Calculate: 


X^U^ul Xtoiuh ficH^U), Pntwl 

Equations (B.S) to {B.8) 




Equation (B.10) 



Calculatei 

Equation {B.11) 



Yes 



y » ^ « 1 



Calculate: 


r„« 273.15 K 
Sc„M. B„(v). p^k). Hsty) 

Equations (B.12) to {B.25) 


No 






Equation (a.26^^^ 




JYes 


Result: 


*CH. JfN2- 'CO- ''CH 



Figure — Flow diagram for computing intermediate data by iteration 
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MchM = - 2,709 328 + 0.021 062 199 //ch(w - 1) ... (B.5) 

^chM^ ^s/[^ch(^ - V (hn^(y)] - [{xh^^h, + -fco^coV^cHf" - 1] - . . (B.6) 

-^Na («) = 1 - -^Ch(«) - JfCOg - -»^H2 " -^CO ' ... (B.7) 



Ai(") = [jfCH(«)A^CH(«) + -^N2{")^N2]ftn,n(v) + (xcx>2A/c^ +J^H!2A^H2 + JfCO^Co)Pm,nM . . . (B.8) 

Table B.I — Values of the constants used in equation (B.I) 

(adjusted to conform with the molar masses and molar calorific Values in ISO 6976:1 975) 





= 285.83 


MJkmohl 


^CO 


« 282,98 


MJkmol-i 


^Hz 


= 28.013 5 


kgkmoM 


A/cOz 


= 44,010 


kgkmol-' 


Mhz 


« 2.015 9 


kgkmohl 


Moo 


= 28.010 


kgkmol-i 


R 


= 0.0831451 


m3barkmo|-i.K-i 




= 22,414 097 


m^kmoH 


Pn(aiO 


= 1,292 923 


kgm-3 



B.2.1 Iteration with the molar heating value Hon (inner loop) 

Equations (B.1) to (8.8) are applied in sequence so as to obtain the first approximation in the uth iteration step. The 
starting values are: 

Wch(« = 0) = 1 000 MJ-kmof-i 

Bn(v = 0) = - 0.065 m3-kmol-i 

The values of the other constants used in equations (B.I ) to (B.8) are given in table B.I . 

The convergence criterion for this inner iteration loop is that the absolute difference between the calculated density 
of the gas at normal conditions pn(u) and the known density (either measured directly or determined from the 
relative density) of the gas at normal conditions pn is less than 10-8, I.e. 

|Ai-A,(«)|<10-6 ...(B.9) 

If this condition is not satisfied, then an improved value of the molar heating value Hch{^)» for use in equations (B.5) 
to (B.B), is calculated using equation (B.IO) as follows: 

Wch(") = Wch(« " 1) + OHcni'^) . . . (B.1 0) 

where 

£>WcHW = [pn-Ph(«)][p(i>)-ft,{u)]'^ ...(B.11) 

Pn(u) being the density value for the current iteration step [comnr^encing with WchC" - "!)!» 

p(D) being the density determined by equations (B.4) to (B.8) using [//ch(" - 1) + 1] as input for the molar 
heating value. 
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When the left-hand side of equation (B.9) is less than this iteration loop is terminated and iteration with the 
second virial coefficient begins. 

B.2.2 Iteration with the second virial coefficient fin (outer loop) 

The intemnediate values jcch(u)« Jcn2(")» 'co and HchM ^rom the preceding heration and the input data xcoz and x^z 
are used to determine an improved value for the second virial coefficient ^n(v) for the whole gas at normal 
conditions. 

The second virial coefficient for the natural gas is given by the following equation: 

'B(T) = jc? + ^^X2B^2 + 2x1X3^^3 + 2x1X4 St4 + SxiXsBtS + 4^22 + ^2^3^ + ^^2^4^24 + 

+ '1^33 + J^|fi44+ 4^55 ...(B.12) 

Some of the terms that are missing in equation (B.12), i.e. 525. B34. etc., have been found not to improve the 
accuracy of the calculation if included and are therefore set at zero. 

Bn(v) = fi(rn) ...(B.13) 
where 

r = rn =273,15 ...(B.14) 

JCi = xch(«) ...(B.15) 

X2^x^{u) ...(B.16) 

xz=xco2 .,.(B.17) 

XA-Xy^ ,..(B.18) 
X5=xco ...(B.19) 

fiii=^'//0+^/rt/^CH(«) + i»//2WcH^M ...(B.20) 
where the coefficients bno^ bfn and b^z are second-degree polynomials as a function of temperature 

B^,^bso(^)'i■b^o(VT^b^o(^)T^'^[bm(0)^bmn^^ 

+ [fr//2(0) + ^//2(1)7' + 6//2(2)r2]/f2^(«) ... (B.21) 

and the second virial coefficients 5^4. fii5» ^221 ^23. B24, ^33. ^34. ^44 and B55 are also second-degree polynomials 
as a function of temperature, tn the general form 

Bij-^b^{0)^bij(fiT^bij{2)T^ ...(B.22) 
The uniike-interaction virial coefficients B\2 and B\z are given by 

fll2 = [072 + ^875 X 1 0-5 (320 - r)2] (ffi ^ + B22)/2 ... (B.23) 

B^z = - 0,865(Bi ,533)^^ ... (B.24) 
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The coefficients in equations (B,21) to (B.24) are given in table B.2. 

Table B^ — Numerical values for the coefficients b[0)^ b{\) znd b{2) in the temperature expansion of the 
second viriai coefficient for pure gases and of the uniike^interactlon vlrial coefficients 
CThe units of B are m^-kmoH when the temperature is in kelvins) 





y 


MO) 


Ml) 




CH 


HO 


"4,254 68x10-1 


-2,865 00x10-3 


-4,620 73x10-6 


CH 


H^ 


8,771 18x10-4 


-5,562 81 xlO-« 


8,81510x10-9 


CH 


H2 


-8,24747x10-7 


4,314 36x10-9 


-6,08319x10-12 


N2 


22 


-1,44600x10-^ 


7,40910x10-4 


-9,119 50x10-7 


CO2 


33 


8.683 40x10-1 


4,037 60x10-3 


-5.165 70x10-6 


H2 


44. 


-1,10596x10-3 


8,133 85x10-5 


-9,872 20x10-6 


CO 


55 


-1.308 20x10-1 


6.025 40x10-4 


-6,443 00x10-7 


CH + N2 


12 


y = 0,72 + 1 .875 x 10-5(320 


-7)2 




CH + CO2 


13 


y = -0,865 






CH + H2 


14 


-5,212 80x10-2 


2,715 70 x10-< 


-2,500 00x10-7 


CH + CO 


15 


-6,87290x10-2 


-2.393 81 X 10-6 


5.18195x10-7 


N2 + CO2 


23 


-3,396 93x10-1 


1,61176x10-3 


-Z044 29X10-6 


N2+H2 


24 


1,200 00x10-2 


0.000 00 


0.000 00 



The value of Bn(v) obtained from equation (B,13) is used to calculate the vth approximation of p^, n using equation 
(B.4). 

Equation (B.6) is then used, in the inverse way to that in which it was used previously, to obtain a value for //s{v), i.e. 

^SM = h(")WcH(" - 1) + f4^4 + X5//5]An.n{v) • • • (B.25) 

where //4(=//h2) an^i //5(=//co) are the molar heating values at 298,15 K of hydrogen and carbon monoxide, 
respectively. The convergence criterion for the outer iteration loop (iteration counter v) is that the absolute difference 
between the measured superior calorific value Hq and the calculated calorific value //s(v) is less than 1 0-4, i.e. 

I^s -//sH<10'^ ...(B.26) 

If this criterion is not satisfied, then the value for Bn(v) determined from equation (B.13) is used as a new input value 
for equation (B.4) and the whole iteration procedure, i.e. the inner iteration loop (iteration counter u), is restarted 
from equation (B.5) using the current values of Hch{^ " 1) and Pm.n(v)- 

When both convergence criterion (B.9) and convergence criterion (B.26) are satisfied simultaneously, the final 
intermediate data for the mole fractions jcch and xf^^ and for the molar heating value Hqh have been determined. 



B.3 Calculation of viriai coefficients 

The second and third viria! coefficients B{7) and C(7) of a natural gas are now determined from the mole fractions 
jrco2 and jch^ (input data) and jtcHi JCN2 and xco (intermediate data) and the molar heating value Hch (see figures 
B,1 and B.3). 
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Input data (ttilrd Iteration loop): 

P. K X CM. X ICQ, Jf Mj. jr 








■0 



Calculate: 



Calculate: 


/-(/♦273.1SK) 








/.jiit«1to3 


1 


Calculate: 









CaUulatei 










equations (838). 




/(A r) 



Equation (B.4D 



Rgure — Flow diagram for compression factor calculation 
B.3.1 Calculation of S(T) 

The second virial coefncient B(1) is calculateci from equation (B.I 2) by the procedure described in B.2.2 for a 
temperature 



r = / + 273,15 



...{B.27) 



B.3.2 Calculation of C(7) 

The third virial coefficient for a natural gas at a temperature T is determined using the following equation: 

C(T) = x^Ciii + ^^\^zCu2 + 3jcfx3Cii3 +3x2x4Cti4 + SjcZjcgCns + Sjc^xfCiga + 6JC1JC2JC3C123 + 

+ 3xix2Ci33 + C222 + 3x2x3C223 + 3x2x|C233 + Jf3<^333 + ^^^0444 . , . (B.28) 



The possible additional terms that are missing in equation (B.28) have been found not to improve the accuracy of 
the calculation if included and have therefore been set at zero. 
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Furthermore, in equation (B.28): 

Cm =CH0+ CHINCH +<^H2«CH . . . (B.29) 

where cho. ^i and ch2 are second-degree polynomials as a function of temperature, viz; 

Cm = CHo('o) + CHo(1)7' + CHo(2)r2 +[cHi(0) + cHi{1)r+ CHi(2)r2];/cH + 

+ [CH2(0) + CH2(1)r + CH2(2)r2]//2„ . . . (B.30) 

as are C^, C333, C444, C1-15, C223 and C2331 viz: 

Cfjk = c^jk(O) + c^trtr + c^i(2)T2 ... (B.31) 

The coefficients in equations (6.30) and (8.31) are given in table B.3. 

Table B.3 — Numerical values of the coefficients c(0), c(1) and c(2) in the temperature expansion of the third 
virlal coefficient for pure gases and of the unlike-interaction virial coefficients 
(The units of C are mS-kmoR when the temperature is in kelvlns) 













c(2) 


CH 


HO 


-3.024 88x10-1 




1.958 61x10-3 


-3.163 02x10-6 


CH 


m 


6,464 22x10-4 




-4.228 76x10-6 


6.881 57x10-9 


CH 


HZ 


-3,328 05x10-7 




2,231 60x10-9 


-3,677 13x10-12 


N2 


222 


7.849 80x10-3 




-3.989 50x10-5 


6.118 70 x10-« 


CO2 


333 


2.051 30x10-1 




3.488 80x10^ 


-8,370 30x10-« 


H2 


444 


1.04711x10-3 




-3,648 87x10-6 


4,670 95x10-9 


CH + CH + N2 


112 


y = 0,92 + 0,001 3 (r- 


270) 






CH + CH + CO2 


113 


y = 0.92 








CH + CH + H2 


114 


y = 1.20 








CH+CH + CO 


115 


-7,367 48x10-3 




-2,765 78x10-5 


3.430 51 X 10-8 


CH + N2 + N2 


122 


y = 0.92 + 0.001 3(r- 


270) 






CH + N2 + CO2 


123 


>-1.10 








CH + CO2 + CO2 


133 


y e 0.92 








N2 + N2 + 


223 


5,520 66x10-3 




-1.686 09x10-5 


1.571 69x10-8 


N2 + CO2 + CO2 


233 


3,587 83x10-3 




8.066 74x10-6 


-3.257 98x10-8 



The other unlike-interaction virial coefficients used are given by 

. ^ijk=yijk\f^iiiCjyCiyi^ ...(B.32) 

where yiji is given by 

yi 12 = >122 = 0.92 + 0.001 3{r - 270) ... (B.33) 

>113 =yi33 =0.92 ...(B.34) 

>i 14= 1.20 ...(B.35) 

>123=%10 ...(B.36) 

Equation (B.32) shows that the tennperature dependence of the unlike-interaction virial coefficients is determined 
essentially by the temperature dependence of the third virial coefficients for the pure components. 
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B.4 Calculation of the compression factor and molar density 

The very last stage in the calculation of the compression factor and the molar density is to solve equations (1) and 
(2) simultaneously for the given value of the pressure p. For the first approximation in the iteration using w, p^^ is 
given by . 

ASnH' = 0) = /?r/p + 5 ...(B.37) 

where the second virial coefficient B is defined by equation (B.12) for a temperature T (see figure B.3). An improved 
value pxrM is then given by 

ASV) = (/?r/4l + BAn(v^-1)+^:/^(H'-1)l ...(8.38) 

where the third virial coefficient C for the mixture is defined by equation (B.28) for a given temperature T. The 
convergence criterion for the iteration using w is that the absolute difference between the calculated pressure p{w) 
given by equation (B.39) and the given pressure p is less than 10-5 [see equation(B.40)]. 

p(iv) = y?rAT,(w)[l + SAn(w) + Cp^(H')] ...(B,39) 
|p-p(>v^<10-^ ...(B.40) 

If this condition is not satisfied, then the current value for the molar density pmW) >s used as the new value prn(v^ 1) 
in equation (B.36) and an improved value of the molar density pmi^) is calculated. 

However, if the left-hand side of equation (B.40) is less than 10-5, the iteration routine is ended, and pm(M »s the 
final molar density prn* "^^^ compression factor is then given by 

Z = 1 + 5A„+CpS ...(B.41) 

NOTE — The mass density can be calculated as follows: 

p = [^/phfair)pZn T^/(p^2t)] ... (B.42) 

Zand being rounded to four places of decimals before being used in the density calculation. 
Report the density to 3 significant figures. 



B.5 Consistency checks on the SGERG«88 method 

The following tests, which provide partial consistency checks on the input data, shall be applied when carrying out 
calculations by the SGERG method. 

a) The input data shall satisfy the following condition: 

d > 0,55 + 0,97xco2 - 0.45JCH2 . • - (B.43) 

b) The intermediate calculated value for the mole fraction of nitrogen shall satisfy the following conditions: 

-COI^sjcn^ ^0,5 ...(B.44) 

J^Na +JfC02 '^O.S ...(B.45) 

c) Furthermore, the internal consistency of the input data for the third iteration loop shall satisfy the condition:- 

d > 0,55 + 0,4x^2 + 0,97;ccO2 " O.^^^Hj ... (B.46) 
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Annex C 

(normative) 

Example calculations 

The following example calculations shall be used for the validation of computer implementations of the SGERG 
method not cited in annex B. The calculations have been carried out using the validated executable programme 
GERG88.EXE. which incorporates the subroutine SGERG.FOR described in annex B. 



Table C.1 — Input data 





Gasi 


Ga82 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


^COz 




0.006 


0.005 


0,015 


0.016 


0,076 


0.011 






0,000 


0.000 


0,000 


0,095 


0,000 


0,000 






0.581 


0.609 


0.650 


0.599 


0.686 


0,644 




40,66 


40.62 


43,53 


34.16 


36.64 


36.58 


Table C.2 — Results (Z-values) 


Cond 


tions 


Gas1 


Gas 2 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


p 

bar 


t 

•c 














60 


"3,15 


0.840 84 


0,833 97 


0,79415 


0.885 69 


0,826 64 


0,854 06 


60 


6,85 


0,862 02 


0.856 15 


0.82210 


0,901 50 


0.850 17 


0,873 88 


60 


16.85 


0,880 07 


0,875 00 


0,845 53 


0,915 07 


0,870 03 


0,890 71 


60 


36.85 


0,908 81 


0,904 91 


0,88223 


0,936 84 


0,901 24 


0,917 36 


60 


56,85 


0.929 96 


0,926 90 


0,908 93 


0.953 02 


0.923 94 


0.936 90 


120 


"3,15 


0.721 46 


0,711 40 


0.643 22 


0,808 43 


0.695 57 


0,749 39 


120 


6.85 


0.759 69 


0,750 79 


0.690 62 


0.83613 


0.738 28 


0.784 73 


120 


16.85 


0.792 57 


0.784 72 


0,731 96 


0.859 99 


0.774 63 


0.814 90 


120 


36,85 


0.844 92 


0,838 77 


0.797 78 


0.898 27 


0.831 66 


0.862 66 


120 


56.85 


0,863 22 


0,878 32 


0,845 54 


0,926 62 


0.872 69 


0,897 49 



These gases are the same as the six gases in IS0 12213-2, annex C, where the complete molar compositions are 
given. 
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Annex D 

(normative) 

Conversion factors 

D.i Reference conditions 

The reference conditions for which the standard GERG-88 virial equation was developed and which the 
SGERG.FOR computer subroutine uses internally are 

Calorific value by combustion at Ti = 298,15 K (r^ & 25 ^'C) 

p« 101,325 kPa 

gas metered at T2 = 273, 1 5 K (tz = 0*0) 
p« 101 ,325 kPa 

the latter set of conditions are also the reference conditions for relative density. 

Considerable care is needed to ensure that correctly referenced inputs are used for calorific value and relative 
density. Several countries normally use the above conditions, but others use alternative conditions. This can easily 
cause confusion, particularly since the unit of calorific-value measurement in each case may still be MJ nrr3. 
Table D.1 is a guide to which of the major intemational gas-trading countries use which reference conditions. 

For those using non-metric units for calorific value (i.e. Btu ft-3), conversion both of units and of reference conditions 
is required. 

The conversion factors used are taken from reference [3]. 



Table D.I — Nationally adopted metric reference conditions for the measurement of calorific value 





ri rc) 


12 rc) 


Australia 


15 


15 


Austria 


25 


0 


Belgium 


25 


0 


Canada. 


15 


15 


Denmark 


25 


0 


France 


0 


0 


Germany 


25 


0 


Ireland 


15 


15 


Italy 


25 


0 


Japan 


0 


0 


Netheriands 


25 


0 


Russia 


25 


0or20 


United Kingdom 


15 


15 


United States of America 


15 


15 


NOTES 

1 In all countries the reference pressure fs 101 .325 kPa (= 1 ,013 25 bar). 

2 t\ is the combustion reference temperature. 

3 12 is the gas-metering reference temperature. 
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D.2 Units and conversion factors for pressure and temperature 

If the input variables p and / are not in the necessary units of bar and *C, then conversions must be made in order to 
use the Fortran implementation. A selection of appropriate conversion factors is given in table D.2. 

Table 0.2 — Conversion factors for pressure and temperature 



Pressure 

;j(bar) = WkPa)]/100 

p(bar)»[p(MPa)]x10 

p (bar) = Watm)] x 1.013 25 

p (bar) = (^psia)]/14,503 8 

p (bar) =x [/^(pssg) -i- 14,695 9]/14.503 8 

Temperature 

'(•C) = nK)- 273.15 
fCC)«(rC'F)-32yi,8 
rCC) = [r("R)]/1,8- 273,15 



D.3 Units and conversion of calorific value and density between reference conditions 

Because both superior calorific value and relative density are functions of the composition of a gas mixture, and 
because the thermophysical properties of the individual components depend upon temperature and pressure in 
individual ways, it is in principle impossible (without knowledge of the composition) to convert the calorific value and 
the relative density, known at one set of reference conditions, to exact corresponding values for any other set of 
reference conditions. 

However, because the relevant reference conditions are always thermodynamically close together, and because 
natural gases do not vary in their composition to any major extent, it is possible in practice to give conversion 
factors which may be applied to any typical natural gas with essentially no loss of accuracy in the converted 
physical properties. 



D.3.1 Conversioh factors for the units 

If the input variable is not in the necessary units MJ nrr3, then a conversion must be made. 

The appropriate conversion factors for calorific values at a given set of reference conditions (see clause D.I) are 
given in table D.3. 

Table D.3 — Conversion factors for calorific value 

Hs (MJ.m-3) « [//s(kWh m-3)l x 3,6 
Ws (MJ.m-3) = [//s(Btu ft-3)l/26,839 2 
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D.3.i2 Conversion between different reference conditions 

If the input variables Hs and d are not at the correct reference conditions (see clause D.I), then a conversion must 
be made. Table D.4 gives the conversion factors appropriate to the reference conditions in table D-1. 

Table 0.4 — Conversion factors for calorific value and relative density 

Calorific value at reference conditions 
ty = 25 •C, f2 = 0 *C, P2 « 1.013 25 bar 

Hs-Hsih^ 0*C,r2= 0 •C, ;?2 = 1.013 25 bar) x 0.997 4 
"s =//s ('t = 15 **C. = 15 **C, ;?2 = 1 .013 25 bar) x 1 ,054 3 
//s = Ws ('i = 60 *F, /2 = 60 *F, p2 = 1 ,01 5 92 bar) X 1 ,053 5 
^ Ws (^1 = 60 /2 = 60 'F, p2 1 ,01 5 60 bar) x 1 ,053 9 

Relative density d at reference conditions 
ra = 0 **C, pz = 1 ,0 1 3 25 ba r 

= (fg = 1 5 •C, p2 = 1 .01 3 25 bar) X 1 .000 2 
rf = d(f2 = 60 *F,p2= 1,015 92 baror 1,015 60 bar) x 1,000 2 



NOTE — p = 1 ,015 60 bar = 14,73 psia (pressure customary in US). 
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Annex E 

(informative) 

Performance over wider ranges of application 

The standard GERG vlrial equation has been comprehensively tested* over the temperature range 263 K to 338 K 
and at pressures up to 12 MPa, with the GERG databank 15] and the Gas Research Institute data t^l for gases within 
the ranges of composition, superior calorific value and relative. density given for pipeline quality gases (see 4.4.1). 
Within these limits, the uncertainties are as given in 4.5. 

Rough estimates of the uncertainties involved in calculations of compression factors for wider ranges of application 
are plotted in figures E.I to E.4 as pressure-composition pbts for nitrogen, carbon dioxide, ethane and propane* 
respectively. 

In figures E.1 to E.4, the performance of the SGERG equation is illustrated up to a maximum pressure of 30 MPa. 
This is for informative purposes only and is not intended to imply a recommendation that the equation be used 
routinely or uncritically above the normal specified limits. The uncertainty limits are dependent upon pressure, 
temperature and composition, and are also strongly affected by the proximity of the phase boundary. The estimated 
uncertainty limits presented below are based upon less comprehensive data, published as a supplement to the 
GERG databank and upon the databanks in references [5] and [9]. 




B ±0.1%rot0.2% 
. C t0.2%toaO.S% 
0 lO.S % to ±3.0 % 

0 I 1 I I I L— 

0 0.1 0.2 0.3 0.4 0.S 

Mole fraction of nitrogen, x^^ 



Figure E,1 — Estimated uncertainty limits tor the calculation of compression factors of natural gases 

with a high nitrogen content 
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the worst-case limits have always been chosen. Dashed lines are used to separate two regbns of estimated 
uncertainty when the experimental evidence is not sufficient to determine the position of the boundary. The detailed 
composition of the gas will have a strong influence on the position of the phase boundary and the user should, 
therefore, make his own phase boundary calculation. 




Rgure ESt — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high carbon dioxide content 
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eiSO 



I 30 



20 



10 



A A7 3*0,1% 

B x0.1%to±0.2% 

C *0.2%totO.S% 

0 *03%tOi3.0% 



S6ERG-8B equation 
r»263Kto 338 K 




o.os 



0.1 



0.1S 0.2 
Mote fraction of ethane. x^2m 



Figure E.3 — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high ethane content 
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30 



20 




10 



A A/st0,l% 

B ±0.1 % ro sO.2 % 

C tO.2 % to *0,5 % 

0 tOS % to ±3.0 % 



S6ERC-S6 ec)uatton 
r « 263 tC to 33a K 



0.05 



0.1 

Hole fraction of propane, x^jhs 



Fig|ure E.4 — Estimated uncertainty limits for the calculation of compression factors of natural gases 

with a high propane content 



The overall results at pressures up to lOMPa and temperatures within the range 263 K to 338 K can be 
summarized as follows. Only gases having mole fractions within the limits given below will have uncertainties within 
± 0,1 %, ± 0,2 % and ± 0,5 %, respectively, within the given pressure and temperature domain. 



Component 


Mole fraction for an uncertainty vrithin 




±0,1 % 


±0,2% 


±0,5% 


Nitrogen 


<0,20 


^0,50 




Cart}on dioxide 


^0.09 


^0,12 


s0,23 


Ethane 


•sO.IO 


^0,11 


«0,12 


Propane 


< 0.035 


«0.04 


< 0,045 
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Annex F 

(informative) 

Subroutine SGERG.FOR in Fortran 

c 

c* ******** ****#**************♦***** ********************************* ****** 

C last update: 07.08.96 M. Jaeschke, J. Sikora 

C 

Q ****«*****«««*«**** GERG-88 VIRIAL EQUATION *********************** 

C SUBROUTINE SGERG 

C 

C . * SGERG' CALCULATES THE COMPRESSION FACTORS OF NATURAL GASES USING 

C A SIMPLIFIED GAS ANALYSIS 

C 

C LEGAL COPIES OF THIS PRC>GRAM MAY ONLY BE OBTAINED FROM THE 

C MEMBERS OF THE GERG WORKING GROUP ON COMPRESSION 

C FACTORS OF NATURAL GAS AS GIVEN IN THE GERG TECHNICAL 

C MONOGRAPH TM 5 (1991). 

C 

C The calculations are based on the following four 

C input parameters for the gas analysis: 

C (Valid rainges:) 

C -1- X3: mole fraction C02 { 0.0 -> 0.3 } 

C -2- HS: calorific value in MJ/ra^3 *) ( 20 -> 48 ) 

C -3- RM: relative density *) ( 0.55 -> 0.9 ) 

C -4- X5r mole fraction H2 ( 0.0 -> 0.1 ) 

C 

C ♦) note : metering at T = 0.0 C , P = 1.01325 bar 

C : combustion at T » 25.0 C 

C 

C Further input parameters used are: 

C • 

C P : pressure in bar { 0 -> 120 ) 

C TC: temperature in degrees Celsius (-23 -> 65 ) . 

C 

C A CALCULATED value for X2« the molar fraction of nitrogen 

C is returned. 

. C 

C The calculated values are: 

C 

C Z : compression factor 

C D : molar density in mol/m**3 

C 

C **** For some compilers the SAVE option has to be set explicitly **** 
C 

C SYNTAX : CALL SGERG (X2 , X3 , HS , RM,X5 , P. TC, 2, D) 
C 

C The coefficients used in this program are conform with the 

C values given in subroutine GAS682, from July 20,1988, appended 

C to report 8807, Van der Waals Laboratory, Amsterdam. 

C 

C J.P.J. Michels St J. A. Schouten 

C August 16, 1991 

C 

C Values for the gas constant, molar masses, calorific values 

C and the density of air conform with ISO 6976 (1995) . 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE SGERG (X2 , X3 , HS , RM , X5 , P , TC , Z , D) 
IMPLICIT DOUBLE PRECISION (A-H,0-2) 

IF( P.LT. 0.0 .OR. P .GT. 120.0) STOP ' PRESSURE OUT OF RANGE* 
IF(TC.LT."23.0 .OR. TCGT .65.0) STOP ' TEMPERATURE OUT OF RANGE- 
CALL SGERG1(P,TC,X2,X3.X5,HS,RM,Z,D) 
RETURN 
END 
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SUBROUTINE SGERGl ( P . TC , Q2 , Q3 , Q5 , QM , RM^ Z , D) 
IMPLICIT DOUBLE PRECISION (A-H.O-Z) 
COMMON /RBLOK/ AMOL,HS 

COMMON /XBLOK/ XI , X2 , X3 , Xll , X12 ,X13 , X22 , X23 . X33 
> ,X5.X7,X15,X17,X25,X55,X77 
COMMON /MBLOK/ GM1R0,GM1R1,GM2,GM3 ,GM5.GM7,FA,FB,RL,T0, H5.H7,R 
HS = QM 
X3 = Q3 
X5 = Q5 

C IF(RM.LT, 0.55. OR. RM.GT. 0.90) STOP 'REL. MASS OUT OF RANGE' 

C IF(X3.LT. 0.0 .OR. X3.GT. 0.30) STOP 'C02 OUT OF RANGE' 

C IF(HS.LT.20.0 .OR, HS.GT.48.6 ) STOP 'CALOR. VALUE OUT OF RANGE' 

IF((0.55+0;97*X3-0.45*X5) .GT.RM)STOP 'CONFLICTING INPUT' 
SM = RM*RL 
X7 = X5*0.0964D0 
X33 = X3*X3 
X55 = X5*X5 
X77 = X7*X7 
BEFF= -0.065DO 
H = 1000. ODO 
AMOL= l.ODO/ (FA+BEFF) 
K = 0 
KK = 0 
1 CALL SMBER(H,SMT1) 

IF(ABS(SM'SMT1) .GT. l.D-6) THEN 
CALL SMBER(H+1.0D0,SMT2) 
DH= <SM-SMT1)/(SMT2-SMT1) 
H = H+DH 
KK = KK+1 

IFCKK.GT.20)STOP ' NO CONVERGENCY #1' 
GO TO 1 
END IF 
Xll = X1*X1 
X12 = X1*X2 
X13 = X1*X3 
X22 = X2*X2 
X23 = X2*X3 
X25 = X2*X5 
X15 = X1*X5 
X17 = X1*X7 
CALL B11BER(T0,H,B11) 
CALL BBER(T0,B11,BEFF) 
AMOL= 1.0D0/{FAi-BEFF) 

HSBER = X1*H*AMOL+(X5*H5+X7*H7)*AMOL 
IF(ABS(HS-HSBER) .GT.l.OD-4) THEN 
K = K+1 

IF(K.GT.20)STOP 'NO CONVERGENCY #2' 
GO TO 1 
END IF 

IF(X2 .LT.-O.Ol .OR. X2.GT.0.5) STOP 'CALC- N2 OUT OF RANGE* 
IF(X2+X3.GT.0.5) STOP •N2 + C02 OUT OF RANGE' 
IF(0.55+-4*X2+0.97*X3-0.45*X5) .GT.RM) 
+ STOP 'CONFLICTING RESULT FOR N2 ' 

Q2 = X2 
T = TC+TO 

CALL B11BER(T,H.B11) 

CALL BBER(T,B11,B) 

CALL CBER(T,H.C) 

CALL ITER(P,T,B.C.V,2) 

D = l.ODO/V 

END 

SUBROUTINE SMBER(H,SM) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /RBLOK/ AM0L,HS 

COMMON /XBLOK/ XI . X2 , X3 . Xll , X12 , X13 , X22 , X23 , X33 
> ,X5,X7,X15.X17,X25,X55,X77 
COMMON /MBLOK/ GMIRO ,GM1R1.GM2 ,GM3.GM5,GM7, FA,FB,RL,T0,H5.H7 .R 
GM1= GM1R0+GM1R1*H 

XI = (HS-(X5*H5+X7*H7)*AMOL) /H/AMOL 
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c 
c 



X2 = 
SM = 
END 



1.0D0-X1-X3-X5-X7 

(Xl*GMl+X2*GM2+X3*GM3-^X5*GM5+X7*GM7)*AMOL 



SUBROUTINE B11BER(T, H, Bll ) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON 7BBL0K/ BR11H0(3), BR11H1(3), BR11H2(3), BR22(3), BR23(3), 
BR33(3), BR15{3), BR17(3), BR55{3), BR77{3), B25 

T2=T*T 

B11=BR11H0{1) + BR11H0{2)*T BRnH0(3)*T2 
+ +(BR11H1(1) + BR11H1(2)*T + BRllHl {3 ) *T2) *H 
+ +(BR11H2(1) + BR11H2(2)*T + BR11H2 (3 } ♦T2) *H*H 

END 

SUBROUTINE BBER(T,B11, BEFF) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /BBLOK/ BR11H0{3), BR11H1{3), BR11H2(3), BR22{3), BR23(3), 
+ BR33(3), BR15{3), BR17(3), BR55{3), BR77{3), B25 

COMMON/ 2ETA/ Z12 , Z13 , Y12 , Y13 , Y123 , Y115 
COMMON /XBLOK/ XI , X2 ,X3 . Xll , X12 ,X13 , X22 , X23 ,X33 



T2=T*T 

B22=BR22(1) -i- 
B23=BR23(1) + 
B33=BR33(1) + 
B15=BR15(1) + 
BSS^BRSSd) + 
B17=BR17(1) + 
B77=BR77(1) + 
BA13= B11*B33 
IF (BA13 .LT. 



,X5,X7,X15,X17,X25,X55.X77 



BR22(2) ♦T 
BR23{2)*T 
BR33(2)*T 
BR15(2)*T 
BR55(2)*T 
BR17(2)*T 
BR77{2) ♦T 

0,0) STOP 



BR22{3)*T2 
BR23{3)*T2 
BR33 (3)*T2 
BR15{3)*T2 
BR55(3)*T2 
BR17 (3)*T2 
BR77(3)*T2 



NO SOLUTION* 
ZZZ=212+{320.0D0-T) **2*1.875D-5 

BEFF = X11*B11 + X12*ZZ2*(B11+B22) + 2 .0D0*X13*Z13*DSQRT {BA13 ) 

> +X22*B22 + 2.0D0*X23*B23 + X33*B33 + X55*B55 

> +2 .0D0*X15*B15 + 2.0D0*X25*B25 + 2 • 0D0*X17 ♦B17 + X77*B77 
END 

SUBROUTINE CBER{T, H, CEFF) 

IMPLICIT DOUBLE PRECISION {A-H,0-Z) 

COMMON /CBLOK/ CR111H0(3), CR111H1(3), CR111H2{3), CR222(3), 
+ CR223(3), CR233(3), CR333(3), CR555{3), CR117(3) 

COMMON /ZETA/ Z12 , Z13 , Y12 , Y13 , Y123 , Y115 
COMMON /XBLOK/ XI , X2 . X3 ,X11 , X12 ,X13 , X22 , X23 , X33 

> ,X5,X7,X15.X17,X25,X55,X77 
T2=T*T 

C111=CR111H0 (1) 

+ 



+ (CRlllHl(l) 
+(CR111H2(1) 



C222 
C223 
C233 
C333 
C555 
C117 



CR111H0(2)*T 
CR111H1{2) *T 
CR111H2{2)*T 
CR222 (2) *T 
CR223 (2) *T 
CR233 (2) *T 
CR333 (2)*T 
CR555(2)*T 
CR117(2)*T 



CR111H0(3) *T2 
CR111H1(3)*T2)*H 
CR111H2(3)*T2)*H*H 
CR222(3)*T2 
CR223{3)*T2 
CR233(3)*T2 
CR333(3)*T2 
CR555(3)*T2 
CR117{3)*T2 



CR222(1) 
CR223(1) 
CR233(1) 
CR333{1) 
CR555(1) 
CR117(X) 
CA112=C111*C111*C222 
CA113=C111*C111*C333 
CA122=C111*C222*C222 
CA123=C111*C222*C333 
CA133=C111*C333*C333 
CA115=C111*C111*C555 
IF (CA112.lt. 0.0 .OR. 
+ .OR. CA123.lt. 0.0 
+ STOP • NO SOLUTION' 
D3REP=1.0D0/3 .ODO 

CEFF=Xl*Xll*Clll +3 .0D0*X11*X2* (CA112 ) **D3REP* (Y12+ (T-270 , ODO) 
+ *0.0013D0) 

+ + 3.0D0*X11*X3 *(CA113)**D3REP ♦YIS 
+ f 3.0D0*X1'X15 * (CA115) **D3REP -YllS 

+ + 3.0D0*X1*X22 * (CA122) **D3REP » (Y12+ (T-270 . ODO )* 0 . 0013D0) 



CA113.lt. 0.0 .OR. 
.OR. CA133.LT. 0-0 



CA122.LT.0, 
.OR. CA115. 



0 

LT. 



0.0) 
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+ +6.0D0*X1*X2*X3*(CA123) **D3REP *Y123 
+ +3.0D0*X1*X33 *(CA133) **D3REP *Y13 

+ +X22*X2*C222 + 3 . 0D0*X22*X3*C223 + 3 . 0D0*X2*X33 *C233 
+ •»-X3*X33*C333 + X5*X55*C555 + 3 . 0D0*X11*X7*C117 

RETURN 

END 



C 
C 



NO CONVERGENCY #3' 



SUBROUTINE ITER{P, T, B, C, V, Z) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /MBLOK/ GM1R0,GM1R1,GM2,GM3,GM5,GM7, FA,FB.RL,T0,H5,H7,R 
RT = R*T 
RTP= RT/P 

V = RTP+B 
KK = 0 

V = RTP*{1-0D0+B/V+C/V**2) 
KK = KK+1 

IF (KK .GT. 20) STOP 
Z = 1.0D0+B/V+C/V**2 
PA = RT/V*2 . 
IF(ABS(PA-P) 
RETURN 
END 

**************** ************** 

BLOCK DATA 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

COMMON /BBLOK/ BR11H0(3), BR11H1(3), BR11H2(3), BR22(3), BR23(3j. 
+ BR33(3), BR15(3), BR17(3), BR55(3), BR77(3), B25 

COMMON /CBLOK/ CR111H0C3). CR111H1(3) , CR111H2(3). CR222(3), 
+ CR223(3), CR233(3), CR333 (3 ) , CR555 (3 ) . CR117(3) - 

COMMON /2ETA/ 212 , Z13 , Y12 , Y13 , Y123 , Y115 

COMMON /MBLOK/ GMIRO . GMIRI ,GM2 ,GM3 .GM5 , GM7 , PA, FB, RL, TO , H5 , H7 . R 
DATA BR11H0/-0.425468D0, 0.286500D-2, -.462073D-5 /; 



,GE. l.D-5)G0 TO 5 



BR11H1/0.877118D-3, -.556281D-5. 
BR11H2/-.824747D-6, 0.431436D-8, 



BR22 / 
BR23 / 
BR33 / 
BR15 /. 
BR17 / 
BR55 / - 
BR77 / - 
B25 / 0 
DATA CRlllHO/ 
CRlllHl/ 



144600D0 
-.339693D0 , 
-.868340D0 , 
-.521280D-1, 
-.687290D-1, 
-.110596D-2, 
-.130820D0 , 
O.O12D0 / 

302488D0 



0.881510D-8 /. 
--608319D-11/, 
-.911950D-6 /, 
161176D-2, -.204429D-5 /, 
403760D-2. -.516570D-5 /, 
-.25 D-6 /, 
0.518195D-6 /, 
0.813385D-4, -.987220D-7 /, 
0.602540D-3, -.644300D-6 /, 



740910D-3. 



.271570D-3, 
,2393810-5, 



0,195861D-2, -.316302D-5 



0.646422D-3, -.422876D-5, 0,688157D-8 



DATA 



CR111H2/ -.332805D-6, 0.223160D-8, -.3677130-11/ 

CR222 / 0.784980O-2, -.398950D-4, 0.611870D-7 / 

CR223 / 0.552066D-2, -.168609D-4, 0. 1571690-7 / 

CR233 / 0-358783D-2, 0.806674D-5, -.3257980-7 / 

CR333 / 0.205130D-2, 0.348880D-4, -.8370300-7 / 

CR555 / 0.104711D-2, -.3648870-5, -467095O-8 / 

CR117 / 0.7367480-2, -.2765780-4, .3430510-7 / 

212 / 0.72D0 /, Z13 / -.865D0/, 

Y12 / 0,9200 /, Y13 / 0.92D0/, Y123 /I.IODO/, 

Y115/ 1.200/ 



DATA 


GM1R0/-2.709328O0/, 


GMIRI/ 


.02106219900/, 


+ 


GM2 


/ 28.0135D0/, 


GM3 


/ 


44.01000/, 


+ 


GM5 


/ 2.0159D0/, 


GM7 


/ 


28.01000/, 


+ 


FA 


/22.414097O0/, 


FB 


/ 


22.71081100/, 


+ 


RL 


/ 1.292923D0/, 


TO 


/ 


273. 1500/ , 


+ 


H5 


/ 285.83D0/, 


H7 


/ 


282.98D0/ 




R 


/O. 083145100/ 









END 
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